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IS.  Abstract 

This  AGARDograph  is  a systematic  collection  of  experimental  results  for  use  of  engineers 
and  designer.,  limited  to  externally  carried,  unpowered,  unguided  stores.  In  the  sections 
dealing  with  store  trajectories,  the  physical  situation  is  described,  the  method  of  solution 
is  indicated,  and  the  final  results  are  presented.  Emphasis  is  on  store  motion  in  traversing 
the  parent  aircraft  flowfleld,  with  criteria  for  safe  separation  and  methods  of  calculating 
flight  path.  Windtunnel  test  techniques  are  described,  and  some  information  on  flight 
testing  given.  Analytical  treatment  is  compared  with  flight  and  tunnel  results.  There  are 
brief  mentions  of  ejection  systems  and  the  effect  of  stores  on  aircraft  performance. 
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SUMMARY 


The  problem  of  achieving  safe,  repeatable  and  predictable  separation  of 
stores  from  aircraft  has  received  considerable  attention  In  recent  years.  There  are  so 
many  possible  combinations  of  store  and  aircraft  configurations,  flight  conditions  and 
release  parameters  that  some  effective  and  reliable  methods  must  be  applied  to  ensure 
favorable  separation  characteristics.  Since  Information  on  the  subject  Is  contained  In 
a large  number  of  Independent  sources,  this  report  seeks  to  collect  the  r alts  In  an 
organized  manner  for  the  convenience  of  engJ  leers  and  designers.  In  order  to  limit  the 
scope  of  the  problem  somewhat,  the  report  is  concerned  only  with  externally  carried, 
unpowered,  ungulded  stores. 

A textbook  approach  would  be  desirable.  In  which  each  aspect  of  the  subject 
is  described  physically  and  formulated  mathematically.  However,  store  trajectories  are 
usually  calculated  by  computer  code.  The  computational  programs  are  lengthy  and  are 
useless  without  voluminous  Instructional  material.  Consequently,  In  the  sections  of 
this  report  dealing  with  store  trajectories,  the  physical  situation  Is  described,  the 
method  of  solution  Is  indicated,  and  the  final  results  are  presented.  Wherever  feasible, 
the  formulas,  equations,  or  numerical  procedures  are  provided,  although  In  moBt 
instances,  a reader  should  obtain  the  original  computer  program  and  users'  manual.  If 
he  Intends  to  carry  out  any  computations. 

The  major  emphasis  in  this  report  is  on  the  motion  of  the  store  as  it 
traverses  the  flow  field  In  the  vicinity  of  the  parent  aircraft.  Criteria  for  safe 
separation  are  presented,  as  well  as  methods  of  calculating  the  flight  path  of  the 
store.  Wind  tunnel  test  techniques  are  described  and  Borne  Information  on  flight 
testing  is  also  indicated.  Furthermore,  each  analytical  procedure  Is  checked  by 
comparison  with  data  from  flight  or  wind  tunnel  test  results. 

Components  of  the  ejection  system.  Including  racks,  ejectors,  fire  control 
systems,  and  connections  are  briefly  mentioned.  References  to  more  detailed  Information 
are  provided.  Similarly,  the  effects  of  stores  on  aircraft  performance,  and  resulting 
struc;ural  and  design  implications  are  only  outlined  here.  Sources  of  more  detailed 
Information  are  Indicated.  Thus  this  report  presents  methods  of  obtaining  trajectory 
Information  and  a practical  guide  to  related  aspects  of  external  store  carriage  and 
separation. 
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LIST  OF  SYMBOLS 


To  Tacllltate  use  of  the  various  papers  cited  in  this  report,  the  notation 
follows  that  of  the  original  authors.  Consequently,  the  symbols  take  on  different 
meanings  in  different  sections  of  the  report.  Where  a symbol  has  multiple  definitions, 
the  appropriate  sections  are  indicated.  Some  symbols  which  are  used  at  only  one  place 
and  defined  in  the  text  nearby  are  not  repeated  in  the  nomenclature  list. 
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body  reference  diameter  (Section  5) 

aircraft  drag  Increment  due  to  etore  installation  (Section  3) 
store  reference  dimension 

diameter  or  thickness  parameter  for  shock  detachment 
force  components  acting  on  store 
gravitational  acceleration  (Sections  4 and  5.3) 
slope  of  source  strength  distribution  (Section  5.5) 
slope  of  source  strength  distribution  at  section  Xf 
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angle  to  body  tangent  (Section  5.1) 

acceleration  parameter 
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source  strength  at  section  x^ 

non-dimensional  source  strength  at  xk 

dynamic  pressure 

flight  dynamic  pressure 

flight  dynamic  pressure  of  store 
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body  radius  at  section  x^ 


4 


R(1,J> 

r,9 


r 

r 

rJ 

rvl 

S 

s 


u».v8.wa 

u,v,w 
u' ,v' ,w' 


“1 


“s»\i»"s 


ur,ue 


WWC 


ul,v,vl,v*"i,v 


>J 

V. 


vr*ve 


s 

vr,v8 

v- 

v' 


V ,v 
x y 


aerodynamic  Influence  coefficient 
cylindrical  coordinates 
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coordinates 

components  at,  wing  control  point  of  flow  field  velocity  perturbations 
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ejection  velocity 
aircraft  flight  velocity 
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store  weight  (Section  A) 
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1.  INTRODUCTION 


FIG.  1.1  STORE  SEPARATION 


During  World  War  I a pilot  or  bombardier  could  simply  toBS  a bomb  safely  clear 
of  his  aircraft,  and  he  might  even  hit  something.  If  he  flew  close  enough  to  It. 
Unfortunately,  this  technique  became  ob.'.olete  with  the  development  of  the  enclosed  cock- 
pit, If  not  before.  The  pilot's  hand  has  been  replaced  by  a parade  of  different  release 
mechanisms,  and  the  simple  small  bomb  hos  been  succeeded  by  an  awesome  array  of  stores. 
Figure  1.2,  for  example,  shows  a typical  attack  airplane  In  Its  loaded  condition;  while 
Fig.  1.3  Illustrates  the  variety  of  stores  that  are  now  being  carried  by  one  airplane. 

All  this  versatility  has  not  been  achieved  without  penalty.  Besides  the 
obvious  degradation  In  aircraft  performance  that  Is  Incurred  by  the  weight  and  drag  of 
the  stores,  other  problems  arise.  If  someone  refers  to  a "store  separation  problem,"  he 
might  be  concerned  with  the  repeatability  of  ejection  cartridge  impulse,  or  the 
structural  integrity  of  the  ejection  mechanism;  or  he  might  be  referring  to  the  Ingestion 
of  rocket  exhaust  gases  Into  aircraft  engine  inlets,  or  to  a whole  range  of  structural 
and  aerodynamic  effects.  In  this  report,  after  a general  description  of  the  store 
ejection  process,  we  will  focus  on  the  following  requirements: 

a.  Compatibility  - The  stores  must  fit  physically,  electrically,  and  aero- 
dynamlcally  onto  the  allotted  spaces  on  the  aircraft. 

b.  Safe  Separation  - When  released,  a store  must  clear  the  airplane  without 
colliding  with  It  or  interfering  with  its  operation,  and  without  colliding  with  another 
store . 


c.  Predictable  Trajectory  - The  stores  must  follow  a repeatable  predictable 
path  from  the  release  point  to  the  target. 

We  shall  attempt  to  assemble  the  latest  information  on  these  subjects. 

Besides  describing  the  problems,  we  will  present  analytical  and  experimental  techniques 
available  for  their  solutions.  The  treatment  will  focus  primarily  on  methods  of 
assuring  the  safe  separation  and  predictable  trajectory  of  aircraft-launched  stores. 

We  will  limit  the  discussion  to  unpowered,  externally  carried  stores.  Further- 
more, we  will  be  concerned  primarily  with  the  motion  of  the  store  relative  to  the  launch 
aircraft.  The  effect  of  the  store  on  the  airplane  performance  will  not  be  considered. 
Neither  the  structural  loads  on  the  configurations  nor  the  aeroelastic  and  structural 
deformations  of  the  airplane,  the  store,  or  the  release  mechanism  will  be  described, 
although  such  considerations  can  be  very  important  from  the  point  of  view  of  the 
structural  integrity  of  the  system.  Aeroelastic  effects  can  result  In  flutter  of  the 
store-airplane  combination,  or  other  undesirable  interactions.  These  broader  aspects  of 
store  separation  phenomena  lie  outside  the  scope  of  the  present  treatment.  Even  elastic 
deformations  which  directly  affect  the  release  conditions  are  not  discussed  in  this 
report.  The  release  process  will  be  assumed  to  result  in  an  initial  position  >,nd  in 


FIG.  l.J  FLUTTER  MODEL 
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initial  linear  and  angular  velocities  of  the  store,  and  we  will  focus  on  its  progress 
as  It  traverses  the  flow  field  In  the  vicinity  of  the  launch  airplane. 

Although  operational  constraints,  Involving  logistics,  maintenance,  handling, 
etc.,  are  of  considerable  Importance  In  the  design  of  store  separation  systems,  these 
aspects  of  the  problem  are  also  oraltteu  from  the  current  paper. 

Several  documents  are  available  In  which  multi-faceted  store  separation  data 
has  been  collectfd.  Aircraft/Stores  Compatibility  Symposia  have  been  held  at  Eglln  Air 
Force  Base,  Florida,  In  1969;  at  Dayton,  Ohio,  In  1971;  and  at  Sacramento,  California, 

In  1973.  The  proceedings  of  these  symposia  (refs.  1-3)  contain  a considerable  body  of 
information  on  all  aspects  of  the  problem.  At  the  present  time  (January  1975)  the  next 
symposium  In  this  series  is  scheduled  for  Arlington,  Virginia,  in  September  1975. 

Another  source  of  information  on  the  aerodynamic  aspects  of  the  problem  la  the 
Proceedings  of  the  AGARD  Conference  on. Aerodynamic  Interference  of  1970  (ref.  9).  The 
aerodynamics  of  store  separation  Is  briefly  surveyed  In  ref.  5 which  contains,  In  addi- 
tion to  a description  of  the  latest  applicable  analytical  and  experimental 
techniques,  an  extensive  biDllography  with  286  entries. 

The  present  paper  begins  with  a brief  description  of  the  mechanical  elements 
used  in  the  carriage  and  separation  of  aircraft-mounted  stores.  Next  we  discuss  some 
aspects  of  compatibility  between  aircraft  and  stores.  In  the  following  sections  on 
"Safe  separation"  and  "Trajectory  prediction,"  methods  are  described  for  calculating  the 
motions  of  the  stores  as  they  pass  through  the  aircraft  flow  field.  Finally,  we  will 
discuss  experimental  verification  by  wind  tunnel  and  free-fllght  techniques. 

2.  THE  STORE  SEPARATION  PROCESS 

2.1  Fire  Control  System 

In  general  terms  a fire  control  system  is  a specialized  computer.  Inputs  are 
target  direction  and  bearing,  the  launch  airplane's  position  and  flight  condition,  the 
aerodynamic  and  Inertial  characteristics  of  the  bomb,  and  Initial  forces  lroarted  by  the 
launcher.  The  computer  digests  this  Information  and  generates  signals  dlre<  .lng  the 
airplane  to  the  proper  position  with  respect  to  the  target  and  Indicating  when  to 
release  the  bomb. 

The  computer  must  continuously  generate  solutions  in  real  time;  but,  since  it 
Is  carried  In  the  aircraft,  It.  Is  of  restricted  weight  and  size.  Therefore,  its  speed 
and  capacity  are  limited,  and  It  cannot  complete  elaborate  finite  difference  flow-field 
calculations.  Hence,  effects  of  the  airplane's  flow  field  on  the  trajectory  of  the  bomb 
can  only  be  represented  In  some  simple  form 

2.2  Pylon-Mounted  Ejection  Equipment 

2.2.1  Pylon-ejector  geometry 

Externally  carried  stores  are  attached  to  the  aircraft  wing  or  fuselage  by  a 
pylon  and  ejector  rack.  The  pylon-ejector  system  must  serve  the  dual  function  of 
supporting  the  store  (securely.  It  Is  hoped)  during  carriage  by  the  aircraft,  and 
releasing  or  ejecting  the  store  at  the  command  of  the  pilot  or  the  fire  control  system. 
Various  ejection  devices  are  described  In  ref.  6.  Figures  2.1  and  2.2  from  that  report 
show  typical  pylons  and  ejectors. 

Major  components  of  the  ejector  rack  Include  the  hooks  which  fit  Into  lugs  on 
the  store  and  hold  It  during  carriage;  the  sway  braces  which  bear  on  the  store  to  give 
It  lateral  support  during  carriage,  the  ejection  mechanism  which  pushes  the  store  away 
at  release,  and  the  electrical  connectors  which  pass  signals  and  power  between  the  air- 
craft and  the  store  or  the  ejection  rack.  These  Individual  components  will  be  discussed 
briefly, 

2.2.2  Electrical  connections 

On  bombs,  mines,  and  rockets,  the  electrical  signals  transmitted  to  the  store 
are  primarily  for  the  purpose  of  setting  the  fuze  and  Initiating  release.  At  the  same 
time,  the  bomb  or  store  may  signal  back  Its  condition  of  readlnes  to  the  airplane. 
Electrical  power  may  be  supplied  by  the  aircraft,  or,  for  some  fui  ctions,  can  be  carried 
in  the  pylon.  It  is  even  possible  to  generate  power  momentarily  Ly  the  physical  nature 
of  the  release  process. 

Naturally,  guided  missiles,  which  will  not  be  considered  in  detail  here,  send 
and  receive  all  kinds  of  guidance,  control  and  readiness  signals  during  their  carriage 
on  the  airplane. 

The  need  to  convey  electrical  signals  and  power  can  create  logistical,  opera- 
tional, and  design  problems.  It  would  be  convenient  to  have  the  same  standardized  con- 
nections available  for  all  stores.  However,  this  constraint  sometimes  Imposes  intolerable 
restrictions  on  the  development  of  new  ordnance,  so  that  standards  are  continually 
changing.  It  then  becomes  operationally  difficult  to  be  sure  to  always  have  the  correct 
electrical  connection  for  each  different  store-aircraft  combination.  Computer-aided 
techniques  for  determining  electrical  interface  requirements  are  described  in  ref.  7. 
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Good  design  can  simplify  the  logistics  requirements  and  at  least  prevent  the  possi- 
bility of  Incorrect  connections. 

Connectors,  themselves,  can  cause  problems.  Under  combat  conditions,  It  Is 
easy  to  make  the  wrong  connection  or  none  at  all.  Therefore  It  Is  desirable  to  make  the 
designs  as  simple  and  foolproof  as  possible.  They  must  also  retain  their  Integrity  In 
spite  of  adverse  conditions  that  might  be  encountered  during  launch  or  carried  flight. 
Some  typical  design  Information  Is  available  In  ref.  8. 

Arming  wires  and  other  pieces  of  connector  hardware  may  break  loose  during 
store  separation.  Such  debris  may  become  entangled  with  aircraft  control  surfaces  or 
Imbedded  In  other  structures,  or  possibly  Ingested  Into  aircraft  engines  with  dangerous 
consequences  (ref.  9). 

2.2.3  Election  or  release  mechanisms 

The  main  functional  component  of  the  support  pylon  is,  of  course,  the  ejection 
mechanism.  One  type  of  ejector  is  illustrated  schematically  In  Pig.  2.3.  It  Is  powered 
by  the  expanding  gases  released  by  an  explosive  cartridge.  The  mechanism  releases  the 
store  by  opening  the  hooks  from  which  It  Is  suspended.  At  the  same  time,  a piston  (the 
design  Illustrated  in  Fig.  2.3  uses  two  pistons)  pushes  on  the  store  giving  it  an 
initial  impulse  to  help  It  clear  the  airplane. 

The  explosive  release  process  subjects  components  of  the  ejector  to  high 
structural  loads  which  can  result  in  failures  of  critical  components.  The  explosive 
gases  also  erode  passages  and  deposit  particles  which  degrade  performance.  Resulting 
design  and  maintenance  problems  are  discussed  In  ref.  10,  for  example.  The  development 
of  a reliable  ejection  system  requires  careful  testing  of  proposed  systems  and  continual 
reporting  and  analysis  of  operational  problems. 

Some  ejectors  have  one  piston  (or  "foot");  others  have  two.  Also,  some  two- 
plston  ejectors  have  Interchangeable  orifices  so  that  the  Impulse  can  be  adjusted 
differently  between  front  and  rear  ejection  feet  to  compensate  for  differences  In  store 
properties  or  flight  conditions.  In  all  current  adjustable  systems,  any  alterations 
must  be  made  while  the  aircraft  Is  on  the  ground. 


Effects  of  different  ejection  force  distributions  are  indicated  In  ref.  11. 
The  explosive  power  la  delivered  by  a cartridge  inserted  Into  the  rack  when  the  stores 
are  mounted  to  It.  Various  cartridge  charge  levels  are  available.  Some  Idea  of  the 
statistical  variation  between  nominally  identical  cartridges  can  be  3een  from  the 
following  data  taken  from  ref.  12. 


average  impulse  (19  tests) 
RMS  variation  in  Impulse 
average  lag  time  (25  tests) 
RMS  variation  in  lag  time 


81.31  lb-sec 
7.8* 

5.12  millisec 
16.7* 


A variation  of  one  millisecond  in  lag  time  (the  time  between  electrical  signal 
and  store  first  motion)  results  in  only  a 1-foot  change  In  range  for  a store  released 
from  an  airplane  flying  level  at  1000  feet  per  second.  However,  If  the  aircraft  Is 
maneuvering,  small  changes  In  lag  time  can  cause  significant  variations  in  store 
Impact  points. 


Any  change  In  total  impulse  will  result  In  a corresponding  change  In  the 
Initial  linear  and  angular  velocities  of  the  store.  Essentially  the  velocity  of  the 
center  of  gravity  of  the  store  is  incremented  as  follows  (neglecting  the  flexibility 
and  finite  mass  of  the  wing  and  support  and  aerodynamic  effects). 


AV  . • I AI  (2.1) 
m 

where  AVej  Is  the  increase  of  store  velocity  at  the  end  of  the  ejection  stroke  due  to  an 
increase  of  total  Impulse  AI  for  a store  of  mass  m.  The  initial  angular  rate,  ue j , will 
Increase  correspondingly 


Aue j - i AI  (2.2) 

where  1 is  the  moment  arm  from  the  store  center  of  gravity  to  the  ejector  foot  location 
and  1 is  the  moment  of  Inertia  of  the  store. 

The  effects  of  such  changes  in  initial  conditions  on  the  subsequent  trajectory 
of  the  store  can  be  assessed  by  making  trajectory  calculations  using  methods  described 
in  Sec.  5. 

2.3  Sources  of  Dispersion 

In  order  to  place  bombs,  or  other  stores  accurately  on  target,  dispersion  must 
be  minimized.  We  will  define  three  sources  of  dispersion  as  follows;  aim  errors,  launch 
errors,  and  ballistic  dispersion.  Ballistic  dispersion  results  from  the  variation  in 
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FIG.  2.3.B  CUTAWAY  VIEW  OF  HARRIER  EJECTOR  RELEASE  UNIT, 

ERU-119 
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flight  path  of  different  samples  of  the  same  store.  In  general,  the  ballistic  dispersion 
of  a well  made  stnble  store  Is  small,  on  the  order  of  2-5  mils.  However,  certain 
factors  can  cause  sizable  Increases.  For  example,  mounting  lugs  and  other  protuberances 
can  cause  the  store  to  trim  and  then  the  trajectory  become!  sensitive  to  the  roll 
orientation  of  the  store.  Flight  te3ts  with  wind  tunnel  aid  analytic  comparisons, 
reported  In  refs.  13-15,  Indicate  the  Importance  of  these  Items.  Mass  asymmetries  could 
have  similar  effects. 

Unstable  flights  were  also  observed  In  these  tests.  Unstable  regimes  can 
result  from  Magnus  effects  or  pitch-roll  resonance.  It  Is  desirable  to  cause  the  3tore 
to  roll  slowly  (by  canting  the  fins)  to  average  out  the  effects  of  asymmetries  caused 
by  manufacturing  variations.  Too  high  a roll  rate  can  cause  Magnus  Instability,  and  a 
low  roll  rate  must  be  carefully  controlled  to  avoid  pitch-roll  resonance.  Stability 
requirements  and  design  Implications  are  described  In  refs.  13-16. 

Some  limit-cycle  types  of  resonance  can  be  excited  only  If  the  store  somehow 
reaches  a high  angle  of  attack.  The  launch  flow-field  Interference  can  excite  nonlinear 
Instabilities  and  result  in  anomalous  trajectories.  Bombs  which  experience  ballistic 
instabilities  can  fall  far  short  of  predicted  impact  points,  and  endanger  friendly  areas 
while  missing  Intended  targets.  Therefore,  smooth  launch  Is  required  to  reduce  ballistic 
dispersion  as  well  as  to  reduce  dispersion  due  to  mal-launch.  Since  a well  designed 
bomb  has  low  ballistic  dispersion,  the  primary  perturbations  of  the  trajectory  o»  cur 
during  Its  passage  through  the  aircraft  flow  field.  This  launch-induced  dispersion  Is 
simply  the  variation  In  Impact  point  thac  results  from  this  portion  of  the  trajectory. 
Even  without  excitation  of  unstable  motion,  launch  disturbances  can  cause  dispersions  of 
from  5-bQ  mils.  Assuming  that  two  Identical  launches  will  result  In  two  Identical 
trajectories,  the  sensitivity  to  small  variations  In  initial  ejection  conditions,  air- 
craft flight  or  loading  conditions  make  predictable,  repeatable  trajectories  difficult 
to  achieve.  At  high  aircraft  speeds,  where  aerodynamic  forces  become  large  In  comparison 
to  Inertial  loads,  the  situation  is  especially  difficult.  Methods  of  reducing  this 
sensitivity  by  adjusting  the  ejection  system  to  produce  a smooth  launch  are  discussed 
in  ref.  17. 


The  final  cause  of  dispersion  is  the  aim  error.  The  sensors  and  fire  control 
systems  which  aim  the  bomb  are  not  going  to  be  discussed  here,  but  it  is  Important  to 
recognize  the  erfect  on  bombing  accuracy  of  tne  factors  relating  to  dispersion. 

Figure  2. A Illustrates  this  point. 

On  this  figure  "ballistic  dispersion”  Includes  both  mal-launch  and  free-fllght 
effects.  If  both  "aim  error"  and  "ballistic  dispersion"  are  small,  then  high  hit 
probability  Is  possible.  Obviously  it  makes  little  sense  to  improve  only  one  component. 
The  numbers  that  would  make  the  plot  quantitative  de>  end  on  warhead  size, target  size, 
range,  and  other  factors.  In  general,  however,  a 5-mll  tolerance  in  both ' trajectory  and 
aim  error  leads  to  a very  effective  system;  while  more  than  10  mll3  in  either  component 
is  not  acceptable  for  accurate  bombing. 


HIT  PROBABILITY 


FIG.  2.4  HIT  PROBABILITY  AS  A FUNCTION  OF 
AIM  ERROR  AND  DISPERSION 


Aim  error  Is  apparently  greater  in 
CJt.bat  th_r.  it.  t„»t  raiige  cehJltltv., 
Improvements  in  aiming  systems  must 
address  realistic  environments. 

The  direct  side  force  control 
system  (ref.  18)  seems  to  promise 
progress  in  that  direction. 

3.  COMPATIBILITY 

3.1  Geometric  and  Functional 

Compatibility 

The  requirement  for 
physical  clearance  is  such  a basic 
one  that  it  may  seem  unlikely  that 
any  possibility  of  interference 
would  Le  overlooked.  However,  an 
airplane  is  an  elusive  configuration, 
.arts,  TuCh  &s  it.nfliiig  gfiW  , fl&pT  , 
control  surfaces,  engine  components 
and  even  wing  tips  are  movable. 

Then  the  geometry  can  deflect 
considerably  under  load,  and  the 
configuration  on  the  ground  supported  by 
the  landing  gear  is  not  the  same  as 
that  in  the  air  supported  by  the 
wings. 

Also  new  stores  and  new 
carriage  and  ejection  equipment  are 
constantly  being  Introduced.  The 
new  geometry  is  not  always  compati- 
ble with  all  of  the  older  components 
with  which  it  might  be  used.  In 
fact,  there  are  so  many  possible 
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combinations  of  stores,  "acks  and  aircraft  that  It  Is  not  su-yrising  that  many  of  them 
are  physically  Incompatible. 

Naturally,  It  Is  desirable  to  check  a new  configuration  on  paper  before  build- 
ing it;  or  to  check  the  fit  on  the  drawing  board  of  a particular  combination  of 
components  before  trying  It  out  with  the  actual  hardware.  Procedures  have  been 
developed  for  performing  this  non-trivial  analysis.  The  one  developed  by  Washmuth 
(refs.  19  and  20)  will  be  described  here. 

The  first  step  is  the  assembly  of  three-view  drawings  showing  aircraft,  racks, 
and  stores  to  the  same  scale.  The  aircraft  drawings  show  cross  sections  at  each  pylon 
giving  interference  lines  with  wheel  well  door3,  flap3,  ailerons,  ground  lines,  etc. 
Examples  of  such  drawings  for  the  Douglas  A-9  aircraft  are  shown  in  Figs.  3.1  and  3-2 
(from  ref.  19).  It  is  sometimes  convenient  to  superimpose  the  interference  lines  of  a 
number  of  aircraft  on  a single  chart.  Then  they  can  be  checked  simultaneously  for 
clearance  problems  with  a new  store.  One  such  drawing  is  shown  in  Fig.  3.3. 

The  next  step  in  the  procedure  in  to  lay  the  Btore  and  rack  drawings  on  top 
of  the  aircraft  layouts  and  identify  possible  interference  in  each  vie-.'.  Such  a super- 
position is  illustrated  in  Fig.  3-9. 

Finally,  the  rack-store  interface  is  checked  with  drawings  that  show  the 
locations  of  ejector  feet,  sway  braces,  electrical  connections,  and  other  Interfacing 
components . 


In  principle,  the  information  could  be  stored  In  computer  programs  which  would 
then  quickly  check  for  interference.  In  practice,  however,  uslrg  the  procedure  Is  very 
simple;  the  hard  part  being  the  acquisition  of  the  required  dravings.  The  Naval  Surface 
Weapons  Center  has  an  extensive  collection  of  these  interference  test  drawings  of  U.S. 
aircraft,  racks,  and  stores.  Many  are  available  in  ref.  20. 

Besides  fitting  geometrically  onto  the  aircraft,  the  store  must  mate  with  the 
systems  that  operate  and  release  it.  Electrical  connections  ha  re  been  alluded  to 
prev?  ly,  and  a computer-aided  design  technique  is  described  ..n  ref.  7.  Other  hard- 
ware l -faces  include  mounting  lugs,  sway  braces  and  ejector  mechanisms.  Some 

requlr  nts  are  established  in  ref.  8. 

Software  considerations  are  also  important.  The  fire  control  system  must 
contain  trajectory  data  for  the  particular  store-aircraft  configuration  which  it  is 
controlling. 

Elaboration  of  these  points  is  not  the  function  of  the  present  paper.  How- 
ever, while  concentrating  on  store  trajectories,  it  is  necessary  to  remember  that  there 
are  other  parts  of  the  system. 

3.2  Aerodynamic  Compatibility 

3.2.1  Aircraft  performance 

A look  at  Fig.  1.2  will  suggest  r.hat  an  airplane's  performance  is  going  to 
suffer  when  it  carries  external  stores.  Several  methods  have  been  developed  for 
calculating  the  drag  penalty  imposed  by  the  stores.  A correlation  procedure  is 
described  by  La  jey  in  ref.  21.  Essentially,  he  expresses  the  incremental  drag  in  the 
form  of  a polynomial  involving  parameters  that  could  affect  the  drag.  Then  the 
coefficients  art  determined  by  fitting  experimental  data. 

For  a single  store,  Lacey  finds  that  lowering  the  store  increases  the  drag; 
interference  effects  are  large  at  M * .9,  but  much  smaller  at  M ■ .8.  The  trends  are 
similar  for  clusters  of  stores. 

Another  correlation  procedure  is  described  by  Berry  in  refs.  22  and  23.  It 
is  based  on  the  more  detailed  report  of  ref.  29.  The  trends  predicted  by  this  analysis 
are  similar  to  those  observed  by  Lacey. 

Gallagher  and  Dyer  (ref.  25)  have  developed  a somewhat  more  comprehensive 
technique,  also  based  on  correlation  of  data.  Their  procedure  has  been  computerized, 
and  covers  subsonic,  transonic,  and  supersonic  speed  ranges.  Some  comparisons  between 
their  predictions  and  experimental  data  are  shown  in  Figs.  3.5  and  3.6 

In  addition  to  the  incremental  drag,  external  stores  add  weight  thus  further 
penalizing  performance.  The  stores  also  will  produce  pitching  moments  which  must  be 
trimmed  out,  again  penalizing  performance.  The  stores  can  be  distributed  in  such  a 
way  that  their  weight  will  not  cause  excessive  shifts  of  the  center  of  gravity  of  the 
airplane . 


The  drag  force  on  the  stores  shown  in  Fig.  1.2  acts  near  the  center  of  gravity. 
The  relatively  high  wing  of  the  A-7  airplane  puts  the  stores  in  a favorable  position 
with  respect  to  the  effect  of  their  drag  on  the  trim  of  the  configuration.  The  drag 
on  stores  placed  under  a low-wing  airplane  will  induce  a moment  requiring  a change  of 
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Store  Installation  loadings: 

(a)  Cut  I • Wing  Pylon-Mounted  Single  Staroa 

(b)  cut  n - Wir«  Pyloo-Mounted  Multiple  Store* 
2)  Data  Involve  tha  following  alrereft-atore 

eetftinetlons  and  flight  conditions: 

OF-1*  Store*! 

a f-ja/b  wc-fa,  KK-ee,  mk-63,  hk-su, 

6 T-S  M-U7,  M-116,  BUK1C/B,  KK-U3, 

2 F-100  KK-57,  MK-61,  LAU-3/A,  AM0-1D 

0 F-lll  ' 

>AJtM  W lancei  0 to  .95 

a-6a 

• A-7  -A-„  Ranee:  26°  to  72.5 


FIG.  3.5  COMPARISON  OF  PREDICTED  VERSUS  ACTUAL  INCREMENTAL  DRAG 


1)  Plotted  4ct«  m for  tlM  fallowing 
Stort  ZnstiUitioQ  1 null  nn  t 

(a)  Can  I - Wing  PyXoo-Nouotod  Singl*  8torM 

(b)  Caae  H - Nine  Pylon -Mounted  Multiple  Store* 

2)  Data  Involve  tha  following  aircraft- ator* 
cceblnatlona  and  fllcht  conditions: 

0 F-4  Store*: 

0 T- 5A/»  / IK-a.  KK-flC,  MR-®*,  «-117r 


m-cu.,  HA-OK,  HU  J./, 

*•118,  Ktt-lC/B,  MB-28,  LAU-3/A,! 
UkU-10/A,  CBU-39,  SU0-13/A, 
Bnllpiv  B,  Walleye,  Phoenix 
Fuel  Tanka  ! 

M Ranee:  0.6  to  2.0  . . 

-A__  Ranee:  20°  to  72.5° 


FIG.  3.6  COMPARISON  OF  PREDICTED  VERSUS  ACTUAL  INCREMENTAL  LIFT 


tr^m.  However,  the  low  wing  configuration  will  generally  put  the  store  trajectory 
further  away  from  the  horizontal  tall  reducing  chances  of  collision. 
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3.2.2  Aircraft  stability 

Besides  Impairing  performance,  external  stores  will  also  affect  the  stability 
and  handling  characteristics  of  the  carrying  aircraft.  Methods  of  computing  inter- 
ference forces  and  moments  are  presented  in  refs.  22-32.  The  procedures  described  in 
refs.  22-26  are  primarily  based  on  correlations  of  e>perimcntal  data.  Such  methods  are 
simple  and  reliable  provided  they  are  applied  to  configurations  which  are  not  outside 
of  the  scope  of  the  data  used  In  the  correlation.  The  other  papers  use  analytical 
methods  based  on  linear  theory.  Wings  and  pylons  are  represented  by  vortex  lattices, 
while  bodies  are  represented  by  source  distributions. 

The  theories  have  been  applied  to  the  calculation  of  the  forces  nn  stores 
attached  to  aircraft.  Figure  3.7  (from  ref.  27)  shows  the  interference  pitching  moment 
on  the  store,  as  a function  of  Its  position,  predicted  by  Fernandes,  compared  with  wind 
tunnel  measurements.  The  store  is  moved  along  its  axis  as  indicated  in  Fig.  3.0,  which 
illustrates  the  configuration.  Comparisons  of  Coble's  analysis  with  measured  normal 
force  and  pitching  moment  on  a wing  mounted  store  are  shown  in  Fig.  3.9  (from  re:.’.  30). 
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Although  Fernandes' calculations  appear  to  show  better  agreement  with  experi- 
ment than  Coble's,  there  is  not  really  much  difference  between  the  methods.  The  results 
should  be  regarded  as  indicating  that  the  theory  is  useful  but  not  infallible, 

3.2.3  Structural  effects 

The  effects  of  store  loading  on  aircraft  structures,  and  the  stresses  on  the 
store  in  carriage  have  been  investigated  for  a number  of  cases  (refs.  33  and  34,  for 
example).  While  conventional  methods  of  structural  analysis  are  applicable  to  the 
calculation  of  the  aircraft  stresses,  the  investigation  of  all  store  loading  and  flight 
conditions  can  be  an  enormous  undertaking.  In  ref.  33,  Brodnax  and  Ripley  describe  a 
"structural  indices"  technique  to  quickly  sort  out  critical  stress  conditions. 


The  problem  of  specifying  store  support  loads  is  addressed  in  ref.  3^  where 
a proposed  Military  Specification  for  store  suspension  equipment  is  analyzed. 


Structural  problems  which  are  particularly  pertinent  to  store-aircraft  compat- 
ibility involve  flutter  criteria  and  effects  of  elasticity  on  the  ejection  process. 
Investigations  of  flutter  have  been  reported  in  refs.  35-38  among  others.  As  with 
structural  loads,  a primary  flutter  problem  is  that  of  investigating  large  numbers  of 
configurations  and  flight  conditions.  Thus,  for  example,  in  ref.  36  Epperson  describes 
a procedure  in  which  flutter  boundary  curves  are  generated  by  a computer  as  functions 
of  stor.  and  flight  parameters.  These  boundaries  are  then  used  to  reduce  the  computa- 
tional labor  required  to  establish  flutter  clearance. 

Feference  38  presents  a different  viewpoint.  Means  of  suppressing  flutter 
are  discussed  rather  than  methods  of  establishing  conditions  for  flutter  clearance. 


In  ref.  39,  Devan  presents  results  of  calculations  of  the  structural  dynamics 
of  store  separation.  Flexibility  of  the  ejection  system  in  one  case  reduces  the 
ejection  velocity  by  about  10  percent,  and  reduces  the  angular  rate  imparted  to  the 
store  by  about  30  percent.  Thus  the  flexibility  effect  can  be  important. 

3.2.4  Aerodynamic  heating 

In  supersonic  carriage,  aerodynamic  heating  can  have  significant  effects. 

IR  domes  are  very  sensitive  to  heating  effects  since  a heated  dome  Increases  the  back- 
ground noise  at  the  IR  detector  and  hence  makes  the  target  more  difficult  to  distinguish. 
Warheads  and  propellants  are  also  heat  sensitive.  Thus  as  carriage  velocity  increases, 
aerodynamic  heating  can  be  expected  to  receive  Increased  attention.  Shock  wave  impinge- 
ment from  aircraft  components  or  reflected  shocks  from  store-generated  disturbances  can 
result  in  regions  of  high  local  heat  transfer  on  both  the  aircraft  and  the  store. 
Consequently  store/aircraft  combinations  present  some  difficult  heat-transfer  problems. 

An  assessment  is  given  in  ref.  0 by  Van  Aken  and  Markarian,  and  some  particular  test 
results  are  described  by  Matthews,  Baker,  and  Key  in  ref.  41. 

3.2.5  Design  considerations 

An  Important  implication  of  investigations  of  aircraft/store  compatibility  is 
that  store-carrying  configurations  should  be  thoroughly  considered  in  the  design  of  the 
aircraft.  In  ref.  42,  two  weapon-configured  designs  are  examined.  One  design  incor- 
porated a MER/TER  weapon  carriage  syEtem,  and  the  other  was  configured  for  "conformal" 
carriage.  The  airplane  geometries  were  considerably  influenced  by  the  store-carrying 
requirements.  For  example,  the  MER/TER  wing  stations  were  efficiently  carried  by  a 
fixed  wing  canard-controlled  airplane;  while  variable  sweep  wings  were  more  desirable 
for  the  fuselage-mounted  conformal  carriage  arrangement.  On  the  other  hand,  the  landing 
gear  was  heavier  in  the  latter  case  because  fuselage  stowage  space  was  not  available. 
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COMPARISON  BETWEEN  PROGRAM  RESULTS  AND  TEST 


Figure  3*10  shows  the  advanced  MER/TER  design  and  the  conformal  carriage  airplane. 

Studies  of  interference  drag  effects  have  indicated  that  "conformal"  carriage 
can  reduce  aircraft  performance  penalties  incurred  by  the  addition  of  external  stores. 
In  conformal  carriage,  the  stores  are  placed  against  the  fuselage,  as  in  the  configura- 
tion pictured  in  Fig.  3.11  (from  ref.  1)3). 

To  test  this  concept,  an  F-4  fighter  airplane  was  modified  for  conformal 
carriage  and  tested  in  a Joint  U.S.  Air  Force-Navy  project.  The  program  is  described 


The  results  of  the  flight  tests,  reported  in  ref.  1)4,  indicated  an  11-percent 
increase  in  range  when  the  empty  rack  configuration  was  replaced  by  a conformal  carriage 
design  with  no  stores.  For  the  two  configurations  loaded  with  twelve  bombs,  the  range 
penalty  of  the  conformal  carriage  configuration  was  18  percent  less  than  for  the 
conventional  loading  arrangement.  Relative  performance  envelopes  are  shown  in  Fig.  3.12 
(from  ref.  1)4).  The  stability  and  handling  of  the  airplane  did  not  suffer  in  the 
conformal  carriage  modification. 
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Subsonic  and  supersonic  bomb  release  tests  were  Included  In  the  program.  As 
reported  In  ref.  A'j,  all  releases  were  successful  with  dispersion  no  worse  than  that  of 
the  conventional  configuration. 

Besides  conformal  carriage,  other  Improvements  In  weapon-carrying  aircraft 
design  may  be  achievable.  For  exam|  in.  stores  ran  carry  lift,  and  hence  might  be  con- 
figured to  reduce  Induced  drag.  Tlr-  entire  aircraft-store  system  could  be  Improved  by 
Integrated  design  teelm  1 quo:-. . 

A.  HAKE  HEI'Ah.V:  i U 

A.l  Th-'  Jal'c  1 '« • [ - 1 1 ■ a T ! r.  Iroblem 

Whenever  a store  Is  released  In  flight,  It  Is  supposed  to  clear  the  carrying 
aircraft  without  hitting  r damaging  it.  In  many  situations,  the  precise  point  at  which 
the  store  Impacts  on  the  ground  Is  not  of  Interest;  the  only  requirement  of  the  separa- 
tion process  Is  that  the  store  does  not  collide  with  the  aircraft.  Such  conditions 
apply  to  Jettison  of  fuel  tanks  or  other  expendable  containers,  or  to  bombing  of  very 
broad  areas. 

Lightweight,  unstable,  low-drag  stores  are  especially  likely  to  present 
separation  problems,  particularly  at  high  dive  angles.  Trouble-prone  flight  configura- 
tions are  summarized  by  H.  Davis  as  follows:  (quoted  In  ref.  A6). 

"Any  store  fired  forward  or  aft  from  airplanes  having  low  set  horizontal 
stabilizers/stabllators  where  the  store  is  fired  within  or  near  the  lateral 
dimension  of  the  stabllizer/stabilator.  . . 

"Any  store  which  exhibits  low  static  stability  or  is  unstable,  particularly 
when  combined  with  a low  moment  of  inertia.  . . 

"Any  store  which  is  marginally  stable  or  unstable  until  fins  pop  out  after 
release 

"Stores  which  change  configuration  immediately  upon  release  (finned  bombs 
released  In  the  retarded  mode)  particularly  if  hardware  Is  released  from  the 
store  at  the  configuration  change.  . . 

"Ejected  stores  In  which  the  ejector  foot  strikes  the  store  at  a position 
remote  from  the  store  center  of  gravity  (empty  rocket  pod,  partially  loaded 
multiple  ejection  racks  and  CUU-AO/AA  flare  dispensers).  The  seriousness  of 
this  situation  is  compounded  by  a low  moment  of  inertia  and  poor  static 
stability.  . . 

"Any  store  ejected  from  the  inboard  shoulder  statlon(s)  of  a multiple  or 
triple  ejection  rack  mounted  on  a pylon  In  close  proximity  to  the  aircraft 
fuselage.  . ." 

A. 2 Criteria  for  Safe  Separation 

A. 2.1  Types  of  criteria 

It  is  unpleasant  to  discover  by  flight  experience  that  a store  collides  with 
the  launch  aircraft  during  separation.  It  is  desirable,  therefore,  to  have  some  means 
of  predicting  safe  separation  boundaries.  In  principle  these  boundaries  can  be  estab- 
lished by  calculating  the  motion  and  trajectory  of  the  store  as  It  passes  through  the 
nonuniform  flow  field  around  the  airplane,  by  making  wind  tunnel  simulations  of  store 
separations,  or  by  determining  formulas  which  will  distinguish  safe  conditions  from 
unsafe  ones. 

Detailed  trajectory  calculations  require  enormous  amounts  of  aircraft  flow 
field  and  store  characteristics  data  because  of  the  large  number  of  possible  combinations 
of  store-aircraft  configurations  and  flight  conditions.  For  the  same  reason  establishing 
safe  separation  boundaries  through  wind  tunnel  tests  is  an  expensive  and  time-consuming 
process.  Therefore  a simpler  procedure  Is  needed  to  establish  approximate  safe  separa- 
tion boundaries.  Borderline  cases  can  then  be  investigated  in  more  detail  if  necessary. 

A simple  reliable  system  developed  by  Covert  (refs.  A6-A9)  will  be  summarized  here. 

A. 2.2  Assumptions 


The  first  assumption  is  that  an  unsafe  separation  will  manifest  itself  in  the 
Initial  motion  following  ejection.  Therefore,  It  is  only  necessary  to  characterize  the 
trajectory  during  a very  short  time  Interval  after  ejection.  It  Is  assumed,  further, 
that  In  this  crlt  cal  Initial  period  aerodynamic  forces  are  constant. 
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STORE  AND  AIRCRAFT  COORDINATE  SYSTEMS 


4.2.3  Coordinate  system 

How  let  u'  take  a look  at  the  coordinate  system  which  Is  shown  In  Fig  4.1. 

The  analysis  Is  carried  out  in  the  x’,  y',  z'  axis  system. 

Two  Important  characteristics  of  the  axis  system  should  be  Identified: 

a)  The  x'-y'  plane  contains  a barrier  (the  pylon)  tnd  hence  a store 
separation  trajectory  13  considered  unsafe  if  any  part  of  the  store  threatens  to  pierce 
this  plane. 


b)  The  coordinate  system  is  fixed  to  the  airplane;  the  x'-z'  plane  being 
parallel  to  the  aircraft  plane  of  symmetry;  the  x'  axis  lies  in  the  plane  of  the  barrier 
(pylon  support).  Since  we  want  to  examine  the  possibility  of  store  interaction  with 
aircraft  components,  it  is  desirable  to  keep  the  axis  system  fixed  with  respect  to  those 
components,  and  hence  fixed  to  the  airplane. 


In  considering  a pylon  Jettison,  the  origin  of  the  axis  system  can  be  moved 
up  to  the  top  of  the  pylon  so  that  the  x'-y’  plane  coincides  with  the  bottom  surface  of 
the  wing  and,  again,  piercing  this  plane  is  assumed  to  be  disastrous. 

4 . 2 . 4 Motion  of  a point 


The  procedure  is  based  on  a consideration  of  the  motion  of  a critical  point  on 
the  store.  Consider,  for  example,  the  vertical  displacement  (in  the  z or  z'  direction) 
of  some  point  on  the  store: 


S z ' 


a°  r 


(4.1) 


where  V0  is  the  initial  velocity  and  a0  the  initial  acceleration  of  the  point.  The 
acceleration  is  assumed  to  be  constant  during  the  short  time,  t,  of  the  separation 
process. 


Separation  will  be  safe  (in  the 
vertical  plane)  if  Sz ^ o for  all  t at  every 
point  on  the  store.  This  condition  will  be 
satisfied  if 

V0  + a0  S.  > 0 (4.2) 

2 

which  will  certainly  follow  if 
V0  > 0 and  a0  > 0 

4.2.5  Vertical  velocity-acceleration  plane 

Thus  we  could  consider  a plane 
containing  all  possible  values  of  V„  as 
ordinates,  and  a0  as  abscissa  (a  velocity- 
acceleration  plane),  and  we  could  identify 
safe  separation  with  that  part  of  the  plane 
in  which  V0  and  a0  are  both  positive. 

(Fig.  4.2) 


FIG.  4.2  VELOCITY-ACCELERATION 
DIAGRAM 
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8 . .’ . 6 ootoj^s  Crl’erlun 

However,  at  the  origin  of  the 
plane,  the  store  has  neither  velocity  nor 
acceleration,  so  It  will  Just  alt  there. 
That  la  not  a cafe  situation.  According 
to  Schoch  (ref.  80),  a store  which  does  not 
fall  a minimum  distance  6a'  In  a critical 
time  tc  la  a potential  source  of  trouble. 
Substituting  this  requirement  In  Eq.  (4,1) 
gives  the  condition 


(4.3) 


Then  the  safe-separation  region  looks  as 
sfnwn  In  Mg.  4.j. 


'overt ' s criterion  assumes  that  a boundary  can  be  established  In  a non- 
dimensional  initial  velocity-initial  acceleration  plane  that  will  distinguish  safe  from 
unsafe  cot.  Jit  lot.;. . This  boundary  Is  to  be  established  experimentally,  but  from  the  above 
discuss  lot.  It  is  apparent  that  In  a non-dimensional  velocity-acceleration  plane,  the 
safe  separation  region  will  be  bounded  by  the  horizontal  axis  and  Schoch's  line,  which 
lr.  nth- dimensional  form  becomes 


where  r Is 
an  i A,  '*r" 


l 

it 


(— )(t0 
£ 


(4.  11) 


maximum  body  radius  and  t0  /—  is  a non-dimensional  critical  time.  At 
ini*  ini  vel  ocities  and  accelerations  lr.  the  vertical  direction  In  'Covert’s 


Figure  4.4  shows  the  safe  separation  region  defined  by  these  criteria  for 
1 j/P*  “ ;*7a,  a Vilue  that  appears  to  agree  with  experiment.  The  circles  on  the  figure 
represent  flight  tests  In  which  the  store  cleared  safely  in  the  vertical  plane. 

tc;Bi  a store  against  this  criterion,  it  Is  necessary  to  show  that  no  point 
on  the  store  pierces  the  safe  separation  boundary,  taking  account  of  the  store's 
angular  motion  as  well  as  its  translation.  Explicit  kinematic  relations  are  given  In 
..ectiop.  4 , 4 , leneral  ly , r.ly  the  nose  and  the  tail  of  the  store  need  to  be  Investigated. 


FIG.  4.4  SAFE  SEPARATION  BOUNDARIES  IN  THE  VERTICAL  PLANE 
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<t.2.7  Safe  lateral  separation 

The  assumptions  underlying  the  vertical  plane  criteria  can  be  applied  to 
determine  analogous  lateral  clearance.  In  a uniform  force  field 

y'*  Ajt  + I i^t2  (4. 5) 

where  and  Ajj  are  Initial  lateral  velocity  and  acceleration,  respectively,  of  some 
x-statlon  along  the  store.  Now  If  the  store  falls  far  enough  to  clear  an  obstacle.  It 
must  travel  the  vertical  distance 


sz<  • ^tc  + | a2tc2 


before  going  the  permissible  lateral  distance  y’.  The  geometry  Is  shown  In  Fig.  4.5. 


FIG.  4.5  LATERAL  COLLISION  BOUNUAKIES 


From  Eq.  (4 . 6) , 


the  time  to  drop  the  distance  &z’  Is 


(*1.7) 


and  then  this  value  Is  substituted  into  Eq.  (4,5)  to  see  If  the  store  will  clear 
(y’<  6y’).  It  should  be  noted  that  If  the  vertical  acceleration  is  negative  (upward), 
then  the  store  will  reach  some  lowest  vertical  position 


z 


max 


at  tc  - 

*2 


If  the  obstacle  extends  below  this  point  (Az’  > zmax) » the  store  will  collide 


with  this  obstacle  unless  the  maximum  lateral  displacement  is  less  than  the  distance  to 
the  obstacle 


A,2 

yfcax  * - «»' 


A 3 and  A 4 


Substituting  tc  from  Eq.  (4.7)  and  letting  y' ■ Ay',  provides  a relation  between 


*3  ■ f^-^c 


(4.8) 


or.  In  non-dimensional  form 


where  t. 


-Aj  tv^Aj2  + 2A  z ' A 2 


A v ' 


3_  . i . 1 

2 6^1 

e 2 


(4.9) 


FIG.  4.6  SAFE  SEPARATION  BOUNDARIES  IN  THE  LATERAL  PLANE 

The  boundaries  for  safe  separation  In  the  lateral  directions  are  shown  In 
Fig.  L . 6 . For  positive  lateral  displacement  (starboard),  the  upper  right  boundary  applies 
and  Sy ' and  is'  (contained  ir.  the  formula  for  t0)  are  the  lateral  and  vertical  clearance 
distances  on  that  side  of  the  missile.  For  negative  displacements,  the  lower  left 
boundary  applies  and  the  normalizing  factors  must  use  6y'  and  Sz’  appropriate  to 
obstructions  on  the  port  side  of  the  missile  (6y'  is  defined  to  be  positive).  The  data 
from  ref.  L7  cannot  readily  be  transferred  to  Fig. '1.0  because  the  normalizations  are 
different. 


To  illustrate  the  effects  of  maneuvers,  consider  the  case  of  an  acceleration 
in  the  vertical  plane.  Assume  that  at  the  instant  of  release  the  airplane  is  flying  in 
-a  circular  path  about  some  center  (Fig.  4.7). 

‘his  maneuver  will  result  in  a relative  displacement  of  a point  on  the  store 
with  respect  to  the  airplane.  The  displacement  will  result  from  a velocity  term  and 
an  acceleration  term.  The  acceleration  term  is  Just  the  centrifugal  acceleration  of 
the  store,  while  the  velocity  comes  about  because  the  airplane  is  rotating  about  its 
own  center  of  gravity,  hence  putting  more  distance  between  itseir  and  the  "store  in 
proportion  to  the  distance  forward  of  Its  center  of  gravity 

The  acceleration  term  is  • 


ng  cos  (80-y0' 
and  the  velocity  is 

x0gg(n-l ) (small) 

where  r.  * number  of  g's  acceleration  which 
the  airplane  is  experiencing. 

Then  the  boundary  (Cchoch's  criterion; 
for  the  maneuvering  case  takes  the  form 


RELATIVE  DISPLACEMENT 


PATH  OF  AIRPLANE 


PATH  OF  STORE 


(~  +N) 


(4. 10) 


FIG.  4.7  MOTION  OF  STORE  RELATIVE  TO 
MANEUVERING  AIRPLANE 
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where  t0  » critical  time,  as  previously  selected  for  Schoch's  criterion;  t0n  « /7?  ■ 

g 

characteristic  time;  N ■ n cos(e0-y0)  - cos  80  cos  $0 


In  a velocity-acceleration  diagram,  on  the  safe  separation  boundary 


is  replaced  by 


eg 


gr 


2 r V- 


(n-1)  where  this  correction  term  Is  generally 


negligible.  In  other  words,  the  line  representing  Schoch's  criterion  moves  more  toward 
the  origin  as  the  pitch  acceleration  Increases. 


*1.2.9  Effect  of  roll  maneuver 


If  the  store  Is  out  on  a wing,  then  it  Is  possible  to  insure  a safe  delivery 
by  rolling  the  airplane  away  from  the  store  Immediately  after  release.  Do  not  roll  the 
wrong  way  II 

4 . 3 Pylon  Jettison 

*4.  3.1  Configuration 

The  Jettison  situation  Is  one  In  which  a safe  separation  criterion  Is  partic- 
ularly pertinent,  because  In  this  case,  once  the  hardware  safely  clears  the  aircraft, 
there  is  little  Interest  in  its  subsequent  trajectory. 

The  Jettisoned  hardware  may  have  any  shape,  but  It  will  be  described  here  as 
a store-pylon  combination. 

*1.3.2  Separation  characteristics 

a)  Such  a configuration  complicates  the  separation  problem  because  It  can 
have  relatively  large  aerodynamic  forces  and  moments  which  give  rise  to  roll,  pitch 
and  yawing  moments. 


b)  The  principal  Inertial  axes  are  Inclined  with  respect  to  the  coordinate 
system,  hence  rolling  motion,  for  example,  can  be  converted  Into  pitching  and  yawing 
motion. 

*1.3.3  Coordinate  geometry 

The  geometry  Is  again  Illustrated  by  Fig.  *4.1.  The  origin  In  this  case  13 
shifted  to  the  top  of  the  pylon  at  its  Junction  with  the  wing. 

The  longer  the  pylon,  the  more  pronounced  will  be  the  unsymmetrlc  effects. 

*t.  3.  **  Analysis 

The  analysis  is  straightforward  but  the  Euler  equations  for  angular  momentum 
lead  to  a set  of  simultaneous  equations.  The  primary  effect  of  the  pylon  Is  to  Induce 
a rolling  moment  and  hence  a roll  rate  ux. 


This  effect  can  be  expressed  In  the  following  form 


rolling  moment  I correction  for  product 


inertia 


of  Inertia  *-crms 


time 


oix 


i(2  _ x 1,2 

l < zz  j K*y 

Hex  k2  -(k4  /k^~ ) 
zz  xy  xx 


. t 


(*4.11) 


where  kxx  and  kzz  are  radii  of  gyration  and  kXy  Is  a product  of  gyration;  1 Is  the 
rolling  moment;  X,  z are  the  coordinates  of  the  center  of  pressure  of  the  pylon. 


This  shows  the  transient  Increase  of  roll  rate  with  time.  Eventually  a 
relatively  steady  roll  rate  will  be  reached,  but  not  '‘•■’ring  the  short  time  for  which 
this  expression  Is  appropriate.  The  roll  will  be  higher  order  in  time  than  other  terms 
and  hence  is  neglected  In  this  short-time  analysis. 


Similarly,  there  Is  a correction  to  the  pitch  plane  motion.  This  correction 
consists  of  an  addition  to  the  pitch  plane  motion  resulting  from  the  yawing  moment 
acting  on  a skewed  principal  axis  of  inertia.  A displacement  in  the  pitch  plane  will 
then  appear. 


This  effect  is  accompanied  by  a corresponding  reduction  of  the  yaw  plane 
motion  due  to  a yawing  moment  since  some  of  the  motion  appears  in  pitch  and  roll. 


*f* 
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4.4  Determination  of  the  Initial  Velocities  and  Accelerations 

4.4.1  dtore  aligned  with  aircraft 

To  test  the  safety  of  release,  the  motion  of  the  store  relative  to  the  airplane  la 
the  quantity  to  be  examined.  To  Illustrate  the  procedure,  consider  first  the  case  In 
which  the  airplane  Is  flying  in  a straight  path  (not  necessarily  a horizontal  one)  with 
the  axis  of  the  store  aligned  with  the  aircraft  velocity  vector.  Then  assume  that  the 
ejector  Imparts  to  the  store,  at  the  end  of  the  ejection  stroke,  an  Initial  linear  veloc- 
ity v'and  an  angular  velocity  u with  respect  to  the  constant  aircraft  flight  path.  The 
ejector  may  push  the  store  In  any  direction;  as  for  example  In  the  case  of  a shoulder 
station  tof  af  multiple  ejection  rack.  The  velocities  may  be  resolved  Into  components  v^, 
v-  and  uy,  uz,  neglecting  any  axial  motion  that  might  have  been  applied.  The  motion  is 
with  respect  to  aircraft-fixed  coordinates,  but  they  are  the  same  as  the  initial  store 
coordinates  In  this  specialized  case. 


The  initial  velocities  at  some  point  x along  the  store  are  then 


Al  " 


4.4.2  General  store  alignment 


- Xuy 

(4.12) 

+ xw'z 

(4.13) 

N'ow  consider  the  more  general  case  in  which  the  store  axis  1b  Initially  displaced 
from  the  aircraft-fixed  x'  direction  by  the  angles  80,  4<0,  and  In  pitch,  yaw,  and 
roll,  respectively.  These  angles  define  the  projections  of  the  missile  axis’  onto  the 
corresponding  planes  perpendicular  to  y',  z',  and  x'.  The  coordinate  system  is  illus- 
trated in  Fig.  4.8 


Assume  that  the  ejection  system  pushes  the  store  away  from  the  airplane  with  an 
initial  velocity  at  the  missile  center  of  gravity  having  x',  y',  z'  components  of  vx', 
vy',  vz',  respectively.  Furthermore,  initial  angular  rates  about  these  axis  will  be 
assumed  to  be  ux',  uy'i  uz'-  At  some  distance  x from  the  center  of  gravity  of  the  store, 
then,  the  Initial  velocity  components  of  a point  on  the  store  axis  become 


A1  • vz'  + cux'  cos  *o  * b“y'  coa  0o 


tan  <l>n 


/ 1+  tan2i|>0+  tan20o 


- xuy' 


/l+  tan2i|<0+  tan2ec 


(4.14) 


l 


FIG.  4.8  MISSILE  AXIS  IN  AIRCRAFT  COORDINATES 


33 


“ vy'  + 8 008  uz'  + 0 8ln  ♦o  “x' 


■ vy'  + xuz'  1 + xwx'  - ta—9^  (4.15) 

• 1+  tan2i(i0+  tan20o  / 1+  tan2<io+  tan26o 

The  Initial  vertical  and  lateral  accelerations  are  obtained  by  resolving  the 
aerodynamic  forces  and  moments  and  the  weight  along  the  z'  and  y'  axes.  The  forces 
divided  by  the  weight  of  the  store  give  the  non-dimensional  accelerations  of  its  center 
of  gravity.  The  moments  divided  by  the  moments  of  Inertial  give  the  angular  accelera- 
tions from  which  the  components  at  a point  x along  the  missile  can  be  calculated.  For 
example.  In  straight  and  level  flight 

r"  008  9°  + [c«*  - ^c.'] 

where  e'  Is  the  pitch  angle  of  the  aircraft  coordinates  with  respect  to  a horizontal 
plane.  0 

The  Initial  linear  and  angular  velocities  Imparted  to  the  store  are  obtained, 
in  principle,  by  applying  the  calibrated  Impulse  of  the  ejector  cartridge  and  the  aero- 
dynamic forces  acting  during  the  ejector  stroke  to  the  Inertia  of  the  store.  However, 
the  impulse  of  the  cartridge  Is  generally  a function  of  the  force  which  It  must  overcome, 
so  that  the  same  cartridge  develops  a greater  total  impulse  if  It  Is  driving  a heavier 
store  or  If  It  is  pushing  the  store  against  a high  in-carrlage  aerodynamic  load.  There- 
fore the  cartridge  must  be  calibrated  over  a range  of  resistances.  Essentially,  the 
force  imparted  by  the  explosive  cartridge  to  an  ejector  piston  is  a function  of  the 
position  and  velocity  of  the  piston.  The  force  is  larger  and  lasts  longer  If  the  piston 
is  moving  more  slowly  as  would  be  the  case  for  a heavy  store. 

Ab  a final  note,  the  flexibility  and  inertia  of  the  ejection  system  may  also 
affect  significantly  the  Initial  velocities  imparted  to  the  Btore.  Devan's  analysis 
(ref.  39)  gives  some  Idea  of  flexibility  effects.  In  Covert's  report  (ref.  49),  he 
points  out  that  a heavy  ejection  system  imparts  a larger  fraction  of  Its  Impulse  to  the 
store  than  does  a lighter  one. 

5.  TRAJECTORY  PREDICTION 

5.1  Determination  of  Aircraft  Flow  Fields  at  Subsonic  Speeds 

5.1.1  Summary  of  methods 

The  trajectory  of  a store  Is  determined  by  the  initial  conditions  and  by  the 

aerodynamic  and  gravitational  forces  which  the  store  encounters.  We  would  like  to  be 
able  to  specify  airplane  geometry,  store  geometry  and  Inertia,  flight  conditions,  and 
ejection  forces,  and  then  have  a computing  machine  tell  us  what  the  trajectory  will  be. 
Such  computer  codes  have  been  formulated  and  are  continually  being  expanded  and  Improved. 
The  first  step  Is  to  describe  the  flow  field  about  the  carrying  aircraft  In  terms  of  the 
airplane  geometry  and  flight  condition.  Next,  the  forces  on  the  store  must  be  determined 
at  each  position  in  the  disturbed  flow  field;  and  finally  the  trajectory  of  the  store 
can  be  calculated  with  these  known  forces. 

The  aircraft  flow  field  can  be  calculated  or  measured  by  probing  the  flow  about 
a wind  tunnel  model.  A complete  set  of  experimental  data  would  require  flow  surveys  at 
a large  number  of  flight  conditions  over  appropriate  ranges  of  aircraft  angle  of  attack 
and  Mach  number.  The  parameter  of  primary  Interest  is  the  local  flow  direction  (upwash 
and  sldewash)  as  a function  of  position  in  the  flow.  Complete  surveys  of  this  type  are 
rare,  but  some  flow  data  Is  available  and  can  be  used  to  test  the  accuracy  of  theoretical 
methods.  In  addition,  there  may  be  buoyant  forces  due  to  pressure  gradients  In  the  flow 
field. 


Se  reral  analytical  techniques  have  been  developed  in  which  the  aircraft  flow 
field  is  represented  In  some  manner,  and  the  passage  of  the  store  through  the  disturbed 
flow  field  Is  then  calculated.  Various  approximate  solutions  are  described,  for  example. 
In  refs.  50-55*  High-speed  computing  machines  make  feasible  detailed  representation  of 
the  aircraft  flow  field  and  step-by-step  numerical  calculations  of  the  resulting  forces 
and  motion  of  a store.  The  analytical  approximations  are  useful  for  rapidly  estimating 
and  comparing  trajectories  and  for  investigating  the  effects  of  various  geometrical  and 
flight  parameters.  In  this  paper,  however,  we  will  examine  in  detail  only  the  -complete 
numerical  procedures  for  store  trajectory  calculation.  Approximate  methods  are  based  o.i 
the  same  principles,  but  apply  simplifying  assumptions  to  reduce  the  computational 
requirements . 


Ooodwln,  Dillenius,  and  Nielsen  (ref.  56)  have  developed  a computational  proce- 
dure for  predicting  atore  separation  trajectories  at  subsonic  speeds,  up  to  the  critical 
speed.  The  first  part  of  the  analytical  process  Is  the  determination  of  the  flow  field 
at  the  store  location.  The  cross-seotlonal  area  of  the  fuselage  of  the  carrying  air- 
craft must  be  specified  as  a function  of  axial  distanoe.  A computer  program  then 
calculates  the  flow  field  for  the  axial  source  distribution  corresponding  to  the  fuselage 
volume.  In  this  computation,  the  fuselage  angle  of  attack  Is  accounted  for  by  the  Beskin 
upwash  associated  with  the  equivalent  circular  cross  sections. 

A later  modification  (ref.  57)  provides  a more  accurate  treatment  for  bodies  of 
non-circular  cross  section.  In  the  modified  calculation,  cross-sectional  shapes  are 
specified  at  a number  of  axial  stations  by  the  locations  of  control  points  on  the  body 
surface.  Then  the  flow  field  due  to  the  sources  which  account  of  the  fuselage  volume  Is 
added  to  the  flow  field  due  to  a selected  number  of  polar  harmonics.  The  coefficients 
of  the  polar  harmonics  are  determined  by  requiring  a least  square  error  In  flow  tangency 
at  the  control  points.  In  this  process,  the  polar  harmonics  account  for  angle  of  attack 
as  well  as  departures  from  circular  cross  section. 

The  modification  also  allows  for  the  introduction  of  engine  inlets.  The  inlet 
cross-sectional  contour  must  be  specified  at  the  selected  axial  stations  along  with  a 
velocity  ratio  which  accounts  for  blockage  effects.  A velocity  ratio  of  one  Implies  that 
the  flow  Is  unaffected  by  the  inlet,  while  a velocity  ratio  of  zero  implies  complete 
blockage,  and  the  Inlet  is  accordingly  treated  as  a solid  body.  For  intermediate  values 
of  velocity  ratio,  the  fuselage  volume  and  the  local  flow  deflections  at  the  inlet 
stations  are  proportionately  adjusted. 

Racks  and  stores  are  handled  in  the  same  manner  as  the  fuselage  In  terms  of 
their  contributions  to  the  aircraft  flow  field. 

The  wing  planform  Is  specified  In  this  computation  by  a number  of  straight  line 
segments.  The  planform  Is  then  subdivided  into  small  area  elements  by  chordwlse  and 
spanwlse  lines.  Specified  angle-of-attack,  twist,  and  camber  distributions  are  then 
applied  to  fix  the  downwash  on  each  area  element  In  the  plane  of  the  wing.  Each  element 
Is  replaced  In  the  computation  by  a horseshoe  vortex.  The  strength.,  of  the  various 
vortices  In  the  lattice  are  determined  by  f 'rcing  the  downwash  due  to  the  entire  vortex 
lattice  to  match  the  downwash  resulting  from  angle-of-attack,  twist,  and  camber  at  control 
points  on  the  area  elements. 

The  wing  thickness  is  represented  by  a distribution  of  sources  In  the  plane  of 
the  wing.  Again  the  wing  planform  is  subdivided  into  a number  of  area  elements  on  which 
the  slope  due  to  the  thickness  Is  assumed  constant.  The  corresponding  source  distribu- 
tion can  be  determined  directly;  and  then  the  velocity  field  induced  by  these  sources 
Is  easily  calculated. 

Pylons  are  treated  exactly  like  wings  except  that  they  are  assumed  to  be  planar 
(no  twist  or  camber)  and  to  have  straight  leading  and  trailing  edges.  The  fuselage, 
stores,  and  racks  are  analyzed  first,  and  the  wing  and  pylons  are  then  analyzed  simul- 
taneously. In  this  way,  flow  angles  Induced  by  the  other  components  can  be  treated  like 
increments  In  wing  camber  and  twist,  and  hence  the  first  order  Interference  effects  are 
accounted  for  In  the  calculation.  In  the  modified  computation  described  In  ref.  57,  an 
Image  of  the  wing  vortex  lattice  is  constructed  Inside  the  equivalent  body  of  revolution 
representing  the  fuselage.  In  this  way  the  wing- fuselage  interference  is  more 
accurately  modeled. 

In  the  Investigation  of  a store  separation  trajectory,  the  flow  field  at  the 
store  location  Is  calculated  In  the  absence  cf  the  store  in  question.  The  accuracy  of 
the  calculatlonal  procedure  Is  Illustrated  by  consideration  of  the  flow  at  the  location 
of  store  #1  In  the  configuration  shown  In  Fig.  5.1.  The  sldewash  and  upwash  distribu- 
tions along  the  store  axis  in  Its  attached  position  and  one  diameter  down  are  shown  In 
Fig.  5.2. 


The  predictions  are  In  reasonab'  •’  good  agreement  with  the  experimental  data. 

The  various  features  are  correctly  represented.  The  accuracy  improves  as  the  store 
moves  away  from  the  airplane  where  the  flow  disturbances  are  smaller. 

At  high  flight  speeds,  the  aero  /ramic  effects  on  the  store  trajectory  become 
Increasingly  more  significant.  Above  the  critical  speed,  shock  waves  appear  which  can 
Impart  strong  perturbations  to  the  forces  on  the  store.  In  addition,  higher  dynamic 
pressures  associated  with  hlghe.  flight  speeds  also  Increase  the  aerodynamic  forces  and 
moments.  Transonic  and  supersonic  flow  field  calculations  are  not  yet  quite  as  fully 
developed  as  the  subsonic  case. 

5.1.2  Flow  field  about  a body  of  revolution 

The  following  sections  will  describe  the  analyses  of  refs.  56-58.  The  terminology, 
nomenclature,  and  results  summarize  their  procedures. 

The  flow  fields  are  determined  for  the  Incompressible  case.  A Prandtl-Qlauert 
transformation  is  then  applied  to  account  for  compressibility  effects.  The  resulting 
relationships  between  the  compressible  coordinates  (x,  y,  z)  and  the  equivalent 
incompressible  space  (x',  y',  z')  are 


FIG.  5.1  LARGE  STORES  WITH  CYLINDRICAL  AFTERBODY  IN  A 
TER  ARRANGEMENT  UNDER  THE  WING 


EFFECT  OF  VERTICAL  POSITION  ON  DOWNWASH  AND  SIDEWASH  UNDER  TER  CONFIGURATION 
(PARENT  AIRCRAFT  - WBPTS2S3) 
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x' 


(5.1) 


y*  ■ y (5.2) 

z'  “ z (5-3) 


The  corresponding  flow  velocities  are  related  by  the  following  expressions 


u 


(5.0) 


v 


V' 

/ 1-Mis 


(5.5) 


and  the  angle  of  attack  relationship  is 


- 


(5.6) 


(5.7) 


These  transformations  are  applied  to  a given  compressible  flow  problem  to  deter- 
mine an  equivalent  li  compressible  flow  configuration.  The  solution  for  the  incompress- 
ible flow  is  then  re-transformed  to  give  the  desired  flow  field  In  the  compressible 
problem.  Thus  It  is  only  necessary  to  determine  solutions  for  the  incompressible  flow 
about  configurations  representative  of  transformed  airplanes  and  stores. 


Figure  5.3  illustrates  the  geometry  of  a body  of  revolution  representing  a store 
body  or  fuselage.  The  sources  are  distributed  In  such  a manner  that  the  boundary  condi- 
tion of  flow  tangency  at  the  body  surface  is  satisfied.  When  the  proper  source  distri- 
bution has  been  established,  the  resulting  non-dimensional  flow  field  velocities  are 
calculated  from  the  following  equations. 


u*(x*,r*) 


N 

l 

k*l 


[<x*-x»)2  + r»2]3/2 

k 


(5.8) 


FIG.  5.3  COORDINATE  SYSTEM  FOR  AN  AXISYMMETRIC  BODY 
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vi(x* ,r* ) - £ 


«k  «•* 


^ [(x#-xj)2  ♦ r*2]3/2 


(5.9) 


The  quantities  have  been  normalized  as  shown  below 


(5.10) 


» . r_  . 

*R  lD 


(5.11) 


v. 


(5.12) 


(5.13) 


«k 


i|!ti2v 
R ” 


(5.1^) 


Here  Ir  Is  the  reference  length  which  is  taken  to  be  the  length  of  the  Incompressible 
body.  vr  Is  the  radial  component  of  velocity,  u is  the  perturbation  velocity  In  the  x- 
dlrectlon.  Qk  la  the  strength  of  the  source  at  the  location  xk. 

The  source  strengths  must  now  be  established.  To  determine  N sources  Qk  where 
k-1,2-  • • N,  the  slopes  of  the  body  are  prescribed  at  N-2  locations;  the  body  is  required 
to  be  closed,  and  the  nose  is  a stagnation  point.  These  conditions  result  in  N linear 
equations  for  the  Qk  as  follows. 

T.nngency  condition  at  the  stations  x*  where  the  body  slope  is  Bj: 


tan  B 


J 


vr<*j,rJ> 

l+u»(xj,r») 


• N-2 


or,  from  (5.8)  and  (5.9) 


tan  B 


J 


2 

k-1  [(xj-x*)2  + rj*2]3/2 


N 

1+  I 


Qk(x.i-Xk) 


J-1,2,'  • 'N-2 


(5.15) 


(5.16) 


k-1  [(xj-xj)2  + rj2]3/2 


The  linear  equations  in  qJ  can  be  made  more  apparent  by  rewriting  (5.16)  in  the  form 


N 


rj-  tan  Bj(xJ-xk) 


tan  Bj  ■ £ Qk 

k‘1  C(xJ-xJJ)2  + r*2]3/2 


J-1.2,- 


■ N-2 


(5.17) 


Closed  body  condition: 


N . 

v':  Qk 

k»l 


(5.18) 


Stagnation  point  condition  at  x*  ■ r*  • 0 


N Qic  ( ~ xk ) 

-i  - : — j— - 

k“1  (xk)3 


«k 


N 
£ 

k*l  x£ 


- 1 


(5.19) 


BETWEEN  ACTUAL  AND  CALCULATED  SHAPES  FOR  TWO  AXISYMMETRIC  BODIES 
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With  these  N equations  the  N values  of  can  be  determined.  Then,  of  course, 
the  Incompressible  flow  field  is  obtained  from  Eqs.  (5.8)  and  (5.9).  The  Prandtl-Olauert 
transformation  finally  gives  the  compressible  flow  field.  Results  of  sample  calculations 
(Fig.  5-1)  show  that  the  body  shape  can,  In  fact,  be  quite  well  matched  by  a suitable 
source  distribution. 

The  angle  of  attack  Induces  upwash  and  sldewash  velocities  on  adjacent  wings, 
pylons  or  other  bodies.  These  velocities  are  calculated  from  the  solution  for  the  two- 
dimensional  flow  about  a circular  cylinder  giving  the  following  induced  velocities  at 
each  x-statlon. 


v(y,z)  ■ — gyZg  g a2W  (5.20) 

(y  +z  ) 

2 2 

w(y,z)  ■-  y ~z a2W  (5-21) 

(y2  + z2)2 

where  a Is  the  local  body  radius  and  W Is  the  velocity  normal  to  the  body  due  to  its 
angle  of  attack: 


W ■ va  sin  a (5.22) 

5.1.3  Non-circular  fuselage 

The  flow  field  about  a non-circular  body  is  given  In  ref.  58.  The  velocity 
potential  for  the  general  case  Is  written  as  the  sum  of  two  parts. 

*c(r,e)  - *e(r)  + *2<r.8>  - S'  (x)  In  r (5.23) 

2ttv«, 


where  Qe  Is  the  three-dimensional  velocity  potential  for  an  equivalent  body  of  revolution, 
that  is,  a circular  body  having  the  same  area  at  each  cross  section  as  the  actual  body. 
The  flow  field  components  for  this  equivalent  body  of  revolution  are  determined  from  the 
solutions  of  Eqs.  (5.8)  and  (5.9).  The  two-dimensional  Inner  potential,  3® 

obtained  by  a superposition  of  polar  harmonics 


MH  an  cosnB  U (x) 

»2(r,e)  ■ t + -2 S'(x)  In  r (5.2>0 

n“l  rn(B)  2iv. 


MH  Is  the  number  of  harmonics  in  the  series.  The  coefficients  are  obtained  by  applying 
flow  tangency  conditions  at  a number  of  angular  positions  equal  to  or  greater  than  MH. 
A least  squares  fit  at  these  points  provides  the  coefficients  an. 


To  carry  out  this  process,  the  body 
cross  section,  assumed  to  have  lateral 
symmetry,  13  divided  Into  equal  angular 
spaclngs  between  6 • 0 and  8 ■ 180  degrees, 
as  shown  In  Fig.  5.5.  The  flow  tangency 
condition  Is  satisfied  by  ♦ ? neglecting 
any  effects  at  the  body  surface  of  the 
other  components  of  the  combined  flow 
potential  (Eq.  (5-23)).  The  crossflow 
components  under  consideration  are  thus 


3<t> 


MH  na  cos  n8 

I — 2 

n»l  rn+* 


Ux(x)S' (x) 

+ 

2ttv«r 


(5.25) 


and 
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MH  na 
Z — 
n-1  r 


n ain  n8 
n+1 


(5.26) 


Neglecting  u,  compared  with  Ux,  the  flow 
tangency  condition  determines  the 
coefficients  an  through  the  relation 


FIG.  5.5  CROSS  SECTION  OF  NON-CIRCULAR  BODY 
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h dv(8) 
v„  dx 


»(x) 
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r ux3'(x) 

L 2»v-r 


MH  nan  cos  n6  *1 

I —r: cos  [g-(<n-  21)] 

n-1  rn+1  J 2 


MH  na_  sin  ne 

♦ r — — — — sin  [e-(m-  1)]  (5.27) 

n-i  rn+1  7 

In  this  equation,  ^ is  the  slope  of  the  body  contour  as  illustrated  in  Pig.  5.6. 
v is  normal  to  the  cross  section  at  r(8),  and  t is  the  corresponding  tangent  vector.  Then, 


- Urn  £ 
djt  Lx*  0 4X 


Ar(8)  cosy 


Lx 


(5.28) 


and  generally  ^ varies  with  8.  m is  the  angle  tn  the  tangent  as  shown.  The  angle  from 
the  vector  r 'to  vector  v is  y and,  as  can  be  Inferred  from  the  figure, 

Y ■ m-8-  \ (5.29) 

and 

m » tan' 

W(x)  is  the  upwash  at  this  cross  section  arising  from  angle  of  attack  or  from  the  flow 
field  due  to  other  components  of  the  configuration.  A sldewash  could  also  be  treated  by 
adding  polar  harmonics  in  sin  n8. 


dr 


£ 

d8 


sin  8 + r cos8 


:os  8 - r sin  8 


(5-30) 


The  coefficients  of  Eq.  (5-27)  are  determined  so  as  to  make  a best  fit  (least 
square  error)  to  the  selected  control  points.  The  number  of  harmonics  and  control  points 
needed  to  represent  the  body  contour  depends  on  the  complexity  of  the  cross  section. 

A fairly  complicated  contour  is  found  to  be  represented  very  well  by  1*1  harmonics  and 
32  control  points. 
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FIG.  5.6 


BODY  CONTOUR  SLOPE 
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Air  Inlets,  at  cross  sections  where  part  of  the  boundary  is  open  to  the  flow, 
are  treated  by  allowing  for  a flow  deflection  through  the  open  portion.  The  tangency 
condition  is  replaced  by  the  slope 


<1- 


v«/dx 


(5.3D 


83  is  the  corrected  slope  of  the  streamline  and  is  directed  toward  the  opening.  vD  is 
the  air  inlet  velocity  produced  by  the  blockage.  The  slope  83  then  substitutes  for  the 

slope  of  the  contour  in  the  calculation  of  the  harmonic  coefficients. 

5.1. 4 Wing  and  pylon  flow  fields 

5 . 1 - ^ . 1 Flow  field  due  to  normal  force 

Flow  fields  about  wing-pylon  combinations  are  generated  by  replacing  the  wing 
and  pylon  by  vortex  lattices  to  account  for  normal  loads,  and  sources  to  account  for 
thickness.  Again,  compressibility  effects  are  accounted  for  through  the  application  of 
a Prandtl-G] auert  transformation  to  obtain  an  equivalent  Incompressible  problem. 

Image  vortices  in  the  fuselage  partially  account  for  uing-body  interference  effects. 

A wc  ng-body-pylon  combination  is  shown  in  Fig.  5.7.  The  fuselage  is  repre- 
sented by  its  equivalent  body  of  revolution.  The  wing  and  pylon  are  represented  by  their 
midplanes  whlc;  have  been  subdivided  into  a number  of  area  elements.  The  wing  may  have 
sweep,  twist,  camber  anc  dihedral  varying  along  the  span.  The  pylon  midsection  is 
assumed  to  be  planar.  Tne  flow  field  due  to  the  wing  and  pylon  normal  force  and  the 
wing  fuselage  interference  is  determined  by  the  velocities  induced  by  a vortex  lattice 
system.  Each  a.-ea  element  is  occupied  by  a horseshoe  vortex  whose  spanwise  segment  lies 
along  the  quarter  chord  of  the  element  and  whose  trailing  legs  lie  along  the  streamwlse 
sides  of  the  element.  Given  the  strength  and  location  of  the  vortex,  the  velocity  field 
which  it  induces  can  easily  be  obtained  from  the  Blot-Savart  law  (see  ref.  59,  for 
example).  Each  horseshoe  vortex  is  accompanied  by  an  image  Inside  the  body  which 

•i 
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FIG.  5.7  WING-PYLON  VORTEX  LATTICE  REPRESENTATION 


FIG.  5.8  HORSESHOE  VORTEX  IMAGING  METHOD  DRAWN  IN  CROSS  SECTION 
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preserves  the  body  contour.  Figure  5.8  shows  a cross  section  of  a general  wing-body 
junction  in  which  the  locations  of  the  tra'.ling  legs  of  the  vortex  image  are  indicated. 
The  origin  of  the  body  coordinate  system  in  at  the  center  of  the  equivalent  body  of 
revolution,  while  the  origin  of  the  wing  coordinates  is  on  the  plane  of  symmetry  at  the 
level  of  the  wing-body  Junction.  The  location  of  an 'image  is  determined  by  the  relation 


r 


iv  ’ 


(5.32) 


where  r„  is  the  radial  distance  to  the  vortex  from  the  center  of  the  body,  and  r<v  is 
the  distance,  along  the  same  radial  line,  to  the  image,  a is  the  body  radius.  This 
combination  of  vortex  and  image,  where  the  strength  of  the  image  is  the  same  as  that  of 
the  external  vortex,  preserves  the  circular  body  cross  section  as  a streamline  in  the 
crossflow  plane. 


Figure  5.9  shows  the  plan  view  of  the  vortex  and  its  image.  The  image  of  each 
leg  starts  at  the  same  streamwlse  (x)  location  as  the  corresponding  external  vortex  leg. 

The  velocity  field  due  to  a general  horseshoe  vortex  is  related  to  the  vortex 
strength,  r,  and  geometry  by  formulas  based  on  the  Biot-Savart  law. 


u(x,y,z)  - k Fu 

(positive  forward) 

(5.33) 

v(x,y , z ) - Fv 

(positive  to  the  right) 

(5.39) 

w(*.y.*>  ■ Fw 

(positive  downward) 

(5.35) 

FIG.  5.9  HORSESHOE  VORTEX  IMAGING  METHOD  SHOWN  IN  PLANFORM 
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The  geometrical  factors  can  be  expressed  as  follows  for  a horseshoe  vortex  on  the  left 
wing  panel. 

Fu  . (z  cos  ♦ - y sin  ♦ ) cos  * 

[x  cos  ♦ - (y  cos  ♦ ♦ z sin  ♦ ) sin  *]2  + (z  cos  ♦ - y sin  #)2 

, | (x  ♦ s tan  »)  sin  » + (y  ♦ a cos  ♦)  cos  » cos  ♦ ♦ (z  ♦ a sin  ♦ ) cos  & sin  ♦ 

* [(x  + s tan  'll)2  + (y  ♦ s cos  $)2  + (z  + s aln  ♦)2]1^2 

_ (x  - S tan  Hi)  sin  it  ♦ (y  - s cos  i)  cos  t cos  t ± (z  - s gin  a)  cos  * sin  ♦ ) 

C ( x - s tan  ili)2  + (y  - s cos  $)2  + (z  - s sin  ♦ )2]1(/2  J 

(5.36) 

p m - z sin  » + x cos  » sin  » 

v ” [x  cos  ♦ - (y  cos  4 + z sin  *)  sin  <i]2  + (z  cos  ♦ - y sin  $)2 


« / (x  ♦ s tan  ♦)  sin  ♦ + (y  + a cos  ♦ ) cos  cos  ♦ » (z  + s sin  ♦ ) cos  i|>  sin  ♦ 
' [(x  + a tan  i| i)2  + (y  + s cos  $)2  + (z+ssln  ♦J2]1^2 


(x  - s tan  ili)  sin  ill  + (y  - s cos  a)  cos  cos  a + (z  - s sin  a)  cos  tfi  sin  a 

222  1/2 
[(x  - a tan  ip)  + (y  - s cos  ♦)  + (z  - s sin  ♦)  ] 
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(z  - s aln  ♦) 

(y  - s cos  $)2  + (z  - s sin  t)2 
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(x  - s tan  ♦) 


[(x  - s tan  a)2  + (y  - s cos  a)2  ♦ (z  - s sin  a)2 


]1/2l 


(z  + s sin  a) 


(y  + a cos  $)2  + (z  + s sin  a)2 
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' [(x  + 


(x  ♦ s tan  ♦) 


[(x  + s tan  ip)2  + (y  + s cos  a)2  ♦ (z  + s sin  ♦J2 


(5-37) 


- x cos  ip  cos  a + y sin  ip 

[x  cos  ip  - (y  cos  ♦ + z sin  a)  sin  <i]2  + (z  cos  ♦ - y sin  ♦)2 

I (x  + s tan  ip)  sin  ♦ + (y  + s cost  ) cos  » cos  a + (z  + s sin  a)  cos  » sin  a 

' [(x  + s tan  ip ) 2 + (y  +scos  a)2  + (z  +ssln  ♦)2]1//2 

(x  - s tan  III)  sin  j ± (v  - a cos  a)  cos  j cos  a + (z  - a sin  ♦)  cos»  sin  a j 

[(x  - s tan  't1)2  + (y  - s cos  ♦J2  + (z  - s sin  ♦ J2]^2  ' 
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(y  - s cos  ♦) 

5 

(y  - a cos  $)  + (z  - s sin 
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(x  - s tan  1(0 

[(x  - s tan  ) 2 + (y  - s cos  $)2  + (z  - s sin  $)2  ]lj^2 


I 
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(y  + s cos  ♦) 

(y  + s cos  $)2  + (z  + s sin  4>)2 


(x  + s tan  i f)  ) 

[(x  + s tan  i|/)2  + (y  + s cos  $)2  * (z  + s sin  $)2]1//2  f 


x,  y,  and  z are  distances  from  the  center  of  the  spanwlse  leg  of  the  horseshoe  vortex, 
as  shown  In  Fig.  5.10.  if  Is  the  sweep  angle  of  this  segment.  In  the  plane  of  the  wing, 
and  ♦ Is  the  dihedral  angle  of  that  part  of  the  wing.  The  sweep  angle  In  the  plan  view, 
i|ip  (see  Fig.  5-10),  Is  related  to  the  sweep  angle  In  the  plane  of  the  wing  by  the 
equation 


FIG.  5.10 


GEOMETRY  OF  TYPICAL  WING  VORTEX  ELEMENT 
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tan  41  ■ tan  4 p cos  ♦ 


(5.39) 


s is  the  semi-span  of  the  vortex  pair  as  shown  on  Pig.  5.10. 

The  same  expressions  hold  for  the  Image  vortices  Inside  the  body,  once  their 
geometric  layout  has  been  established  In  accordance  with  Eq.  (5.32).  Similarly,  the 
pylon  contributions  are  determined  in  the  same  way  with  4 ■ go  degrees.  It  Is  necessary, 
however,  to  first  determine  the  strengths  of  the  vortices. 


Vortex  strengths  are  calculated  by  requiring  no  flow  through  the  wing  at  a 
number  of  control  points.  The  control  points  are  located  along  the  center  span  of  each 
area  element,  three  quarters  of  the  distance  from  the  leading  edge  of  the  element.  The 
vortex  goes  through  the  quarter  chord  and  the  control  point  Is  at  the  three-quarter  station 
because  these  conditions  give  the  correct  lift  and  pitching  moment  on  a two-dimenBional 
airfoil. 

At  each  control  point,  the  normal  velocity  due  to  angle  of  attack,  wing  twist 
and  camber,  and  the  velocity  Induced  by  other  components  such  as  stores  or  fuselage,  is 
cancelled  by  the  normal  velocity  induced  by  all  of  the  wing  and  pylon  vortices  and  the 
body  Image  vortices.  For  each  of  M area  elements  on  the  wing,  these  boundary  conditions 
are  given  by 
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wi.  v 

+ * ■)  cos  * 


(5.40) 


v • 1,2,-  • -,K 

For  each  of  MP  control  points  on  the  pylon,  the  boundary  conditions  can  be  written 

M r M+MP  . V , , , 

- : _n_  (Fv  + Flv  ) - Z -^2-  F - -ixV.  (5.41) 

n»l  4nv.  v.n  ±Yvin  n-M+l  tnv„  Vv.n 

The  symbols  F„  and  Fw  denote  the  geometrical  factors  from  Eqs.  (5.37)  and 
Vv,n  wv,n 

(5.38)  for  the  Influence  of  the  horseshoe  vortex  n at  the  control  point  v.  FlwVin  and 

Fi„  are  the  functions  for  the  corresponding  Image  vortices,  also  obtainable  from 

lvv,n 

Eqs.  (5.37)  and  (5.38).  a Is  the  angle  of  attack  of  the  configuration  and  is  the 

angle  due  to  wing  twist  and  camber  at  control  point  v.  $v  is  the  dihedral  angle  at  the 

control  point.  uljV,  v1>v,  and  w1  v,  are  the  components  of  velocity  perturbation  at  the 

control  point  induced  by  wing  and  pylon  thickness  and  from  other  components  such  as  fuse- 
lage, rack,  and  stores.  Solution  of  these  M+MP  relations  determines  the  values  of  the 
vortex  strengths,  rn.  The  velocity  field  can  then  be  determined  through  Eqs.  (5.36), 
(5.37),  and  (5.38). 

5.1.l(. 2 Wing  and  pylon  thickness 

The  velocity  potential  increment, induced  by  an  array  of  sources  of  uniform 
strength  is  given  by 


v_ 


1_  (\f  Ya  tan 

/ (XrX)2  + 


tan  dXa  dYx 


(5.42) 


(Y-Y1)2  + (ZrZ)2 


kWk 


FIG.  5.11  COORDINATE  SYSTEM  FOR  WING  THICKNESS  SOLRCE  STRIPS 


The  sources  are  distributed  at  all  locations  Xj,  Yj,  Z.,  lying  In  the  area  bounded  by 
the  region  shown  In  Fig.  5. 11.  For  a very  narrow  swept  spanwlse  strip  In  a fixed  vertical 
plane  Zlf  a simplifying  approximation  can  be  applied 


Xx  - X0  - Yi  tan  (5. ^3) 


where  Xj  Is  the  sweep  angle  of  the  strip.  Then  the  integration  results  in 

0AX.  . . 1/2  ■' 

’ TT  008  *1{  ln  [ <CX°  + Va  tan  *1  ' X]Z  + [Y  ‘ Ya]2  + ^Z1  - > 

y I 

- eosa|(|  + (X  - X0)  sin  + Y cos  - tn  [ (tXc  + Yb  tan  *1  - X]2  + [Zj-Z]2  ) 1/2  1 


cos 


+ (X  - X0)sln  + Y cos  **]} 


(5.44) 


AXi  Is  the  width  of  the  strip  ln  the  chordwlse  direction;  Ys  and  Yb  are  the  spanwlse 
limits  (negative  on  the  lift  wing  panel)  shown  ln  Fig.  5.11,  from  ref.  56,  4m  is  the 


J 


sweepback  angle  of  the  strip;  and  6 la  the  slope  of  the  wlrg  thickness  at  the  middle 
of  the  strip. 
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The  perturbation  velocity  components  can  be  obtained  by  differentiating  the 
velocity  potential.  The  resulting  perturbation  velocities  due  to  one  strip  are  then 

4u  . 3 ,44.  64Xi 

^ n W ' IT  008  *i 


( -(Xc  - X ♦ Ya  tan  »x)  + B sin  -(Xe  - X + Yb  tan  + D sin  ) 

| AB  “ CD  ) 

B cos  Y - Yb  + D cos  ^ } 


4 v 
V. 


a aa  04X1  (Y  - Ya  + 


CD 


£ ■ fz  <&>-  ^ cos  | 


-a1  - z>  -(zx  - z) 


AB 


CD 


(5.115) 


where 


/(Xc  - X + Ya  tan  i^)2  + (Y  - Y )2  + (Z!  - Z)2 


COS  ^ 


+ (X  - Xc)  sin  ^ + Y cos 


B - / (Xc  - X Ya  tan  n)2  + (Y  - Y&)2  + (Zj.  - Z)2 


i 

1 

C » / (Xc  - X + Yb  tan  i^)2  + (Y  - Yb)2+(Z!  - Z)2 


*b 

" cos  ^ + (x  - xc>  Bln  *1  + Y 000  *1 


V - / (X0  - X + Yb  tan  *i)2  + (Y  - Yb)2  + (Zj.  - Z)2  (5.40 


Similar  expressions  can  be  derived  for  the  corresponding  source  strip  on  the  right  wing 
panel.  The  resulting  expressions  are  Identical  to  Eqs.  (5.45)  and  (5.46)  except  that 
Ya  and  Yb  are  Interchanged,  since  Ya  is  still  the  Inboard  side  of  the  strip,  and  tan  i|ji 
is  replaced  by  -tan  and  cos  '('i  Is  replaced  by  -cos  4>i.  On  the  right  wing  panel 

Ya  and  Yb  are  positive.  At  a field  point  (X,Y,Z)  the  perturbation  velocity  components 

1A — , !i-  Induced  by  wing  thickness  are  then  given  by  Eqs.  (5.45)  and  (5.46)  summed 

V«  Voo  Voo 

over  all  wing  source  strips  on  both  the  left  wing  panel  t id  right  wing  panel. 
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FIG.  5.12  COORDINATE  SYSTEM  FOR  PYLON  THICKNESS  SOURCE  STRIPS 


The  effects  of  pylon  thickness  are  obtained  in  a similar  manner.  The  pylon  is 
at  a spanwise  location  Yp.  The  geometric  arrangement  is  shown  in  Fig.  5.12,  from  ref. 
56,  and  the  resulting  velocity  increments  are  given  below. 


Au 


cos  l|l^ 


J-XPg  + X,+  Zft,  tan  ^ + F sin  -XPt!  U»Zh  tan  »i  + H sin  »i  1 
' | EF  " QH  ) 


^ cos 

v.  2k  1 ( EF  OH  | 


Aw  9iXl  , ( -za  + 

v-  ■ — 008  ♦ij  


Z - F cos  -Zb  + Z - H cos 


EF 


OH 


(5.47) 
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where 


- A 


(Xpe  - X - za  tan  *i)2  + (Yp  - Y)2  + (Za  - Z)2  + 


COS  Pj 


- (Xp^  - X)  sin  ip^  - Z cos  ip^ 


F ■ /(XPc  - X - Za  tan  <Pi)Z  + (Yp  - Y)2  + (Za  - Z)2 


Q 


X - Zfc  tan  *i)2  + (Yp  - Y>2  + (Zb  - Z)2  + 


2b 

cos  ipj 


- (Xpc  - X)  sin  - Z cos  ^ 


H - /(Xpo  - X - Zb  tan  P^)2  + (Yp  - Y)2  ♦ (Zfc  - Z)2 


(5.18) 


5.1.5  Store  and  ejector  rack  flow  models 

Stores  and  ejector  rack3  are  replaced  by  equivalent  bodies  of  revolution.  The 
flow  fields  about  these  components  are  then  determined  by  the  same  formulas  as  previously 
given  for  circular  fuselages. 

5.2  Force  and  Moment  Calculations  at  Subsonic  Speeds 

5.2.1  Forces  on  store  body 

The  next  step  In  the  process  of  refs.  56-58  for  determining  the  trajectory  of  a 
store  is  to  calculate  the  forces  on  it  at  each  position  of  Its  flight  path.  The  coor- 
dinate system  for  the  force  calculation  is  shown  In  Fig.  5.13.  Delenlus,  Goodwin  and 
Nielsen  calculate  the  forces  In  the  xB,  ys,  zs  coordinate  system  fixed  to  the  missile 
nose.  The  flow  field,  however,  has  been  obtained  In  an  aircraft-fixed  coordinate  system; 
5,n,C.  The  relation  between  the  5,n,C  system  and  corresponding  missile  coordinates 
x,  y,  z is  Illustrated  in  Fig.  5- 11.  The  missile  coordinates  are  similar  to  the  force 
coordinates  except  that  x points  forward  while  xs  points  aft,  and  z Is  positive  downward, 
while  zs  is  positive  upward. 

As  shown  In  Fig.  5. 11,  the  5,n,t  axes  can  be  aligned  with  the  x,  y,  z axes  by 
first  rotating  through  the  yaw  angle  Y about  the  5 axis.  Then  the  system  Is  pitched 
through  the  angle  Q about  the  newly  located  q axis.  Finally  the  system  is  rolled  through 
the  angle  p about  the  x axis.  These  rotations  lead  to  the  relations. 


5 ■ x cos  6 cos  ip  + y(sin  p sin  6 cos  ip  - cos  p sin  p)  + 


z (cos  P sin  0 cos  ip  + sin  p sin  ip) 


n ■ x cos  8 sin  ip  + y(sin  P sin  0 sin  ip  + cos  p cos  tp)  + 


z(cos  p sin  8 sin  ip  - sin  p cos  ip) 


£ ■ -x  sin  8 + y sin  p cos  8 + z cos  p cos  0) 


(5.19) 


The  components  of  velocity  in  tha  store  coordinate  system  are  determined  by 
summing  the  contributions  from  the  free  stream,  the  perturbations  due  to  the  aircraft 
components,  and  the  damping  velocities  induced  by  the  angular  motion  of  the  store.  Thus 
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u„  * U-.  . + u- 

*8 


v-8,ya  + va  + r(x»  ■ x«.m) 


be 


s.z. 


> wg  + q(xg  - x8>m) 


(5-50) 


The  components  of  the  free  stream  as  seen  in  the  store  coordinate  system  will 


U«  ■ cos  6 cos  if  (V„  cos  af  + £)  + n cos  6 sin  * - sin  6 (V„  sin  ar  + c) 
s»xs 

V_.  ■ -(sin  if  sin  6 cos  if  - cos  $ sin  *)(V_  cos  ar  +C) 

8IJS 

-q  (sin  $ sin  0 sin  ip  + cos  * cos  *) 

- sin  * cos  0 (V„  sin  af  + c) 

W«  ■ (V„  cos  af  + 5)(cos  * sin  0 cos  * + sin  * sin  *) 

s » z3 

+ n ( cos  * sin  0 sin  if  - sin  * cos  *) 

+ <V„  sin  af  + c)  cos  * cos  0 (5. 51) 

The  notation,  U„0  _ , denotes  the  component  of  free  stream  velocity  in  the  x.  direction 
s,xa 

of  the  store  coordinate  system.  V„  and  W„  have  analogous  connotations.  af  is 

8,ys  s,ys 

the  angle  of  attack  of  the  parent  aircraft  which  1b  flying  at  the  uniform  speed  V„. 
i,  n,  and  t are  the  components  of  the  angular  rates  of  the  missile  in  the  aircraft- 
oriented  inertial  system. 

The  terms  us,  v8,  and  ws  in  Eq.  (5.50)  are  the  velocity  perturbations  induced 
by  the  fuselage,  wing,  pylon,  rack,  and  other  stores  as  calculated  by  the  previously 
described  methods.  At  each  point  in  the  trajectory  all  of  the  velocity  perturbations 
are  summed  in  the  fuselage  coordinate  system  of  the  incompressible  transformed  space 

giving  u£,  v^,  and  w^  in  the  5,  n,  C coordinates.  When  transformed  back  into  the 

compressible  space,  these  components  become 

Ur  • -4- 
4 B‘ 


w' 

"C  * f 

Then  finally,  in  the  missile  coordinates 

ug  ■ -u^  cos  0 cos  ♦ - vn  cos  0 sin  * + w?  sin  0 

vg  ■ u^(sin  * sin  0 cos  if  - cos  $ sin  if) 

+ vn  (sin  if  sin  0 sin  * + cos*  cos  *)  + w^  sin  * cos 


(5.52) 
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wg  • -ug(cos  ♦ Bin  8 cos  4>  + sin  ♦ sin  ♦) 

-vn(cos  ♦ sin  8 sin  ♦ - sin  ♦ cos  ♦)  - w^  cos  4 cos  8 (5.53) 

The  tern  r(x#  - x.  m)  in  Eq.  (5.50)  accounts  for  the  damping  due  to  yaw,  and 

the  term  q(x8  - x8(m)  is  tn4  pitch  damping  term.  The  symbols  are  defined  in  Pig.  5.14. 

In  the  force  and  moment  calculations,  the  velocities  are  non-dlmenalonalized 
by  the  free-stream  velocity  of  the  store.  Thus 


where 


V„B  ■ [(V«  cos  of  + C)2  + n2  ♦ (V„  sin  af  c )2]1/2  (5.55) 


Body  forces  and  momenta  are  calculated  on  the  basis  of  slender  body  theory. 
The  resulting  expressions  are 


<Vsb  * ^jr*8'0^  <a2w;)dxs  (5-56) 


w 

£ 


(5.54) 


(Cy)SQ  » 


d 

dxs 


(a2v*)dxg 


(5.57) 


(cm) 


SB 


s,o  (x 


s,m 


xs>  ~ (a2“S)d*s 


(5.58) 


(°n)sB 


r8>0(xs,m  - xs>  sf-  U»v»)dx8 


R R 


(5.59) 


Sr  Is  the  reference  area,  and  lR  is  the  reference  length  on  which  the  coefficients  are 
based.  The  upper  limit  of  the  integration,  xs  0 is  the  station  at  which  vortex  separa- 
tion of  the  crossflow  about  the  body  begins,  from  a correlation  by  Hopkins  (ref.  60), 
this  point  can  be  chosen  as  follows 


* -378  + .527  (5-60) 

s 3 

Here  4.  is  the  body  length,  and  xBji  is  the  location  on  the  body  of  the  maximum  negative 
rate  or  change  (boattail)  of  cross-sectional  area. 

In  the  region  of  vortex  separation  a simple  crossflow  theory  is  used  to  provide 
the  body  forces  and  moments. 


»K 


dx. 


(5.61) 


s,o 


(CY)  ■ °dg  f 8 avjvjjdx 
CF  SR  •/xs,o  8 


(5.62) 


^cm^cF  " SRtp  fx  avcws  ^xs,m  “ xs^  dxs 


(5.63) 


S,0 
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where 


f Jt  g 

<cn>CF  ‘ iTj  avMS  <xs,m  - xs}  dxs 

R R xs,o 


v*  ■ / (vj)2  ♦ (wj)2 


( 5 . 6^4 ) 


(5.65) 


and  cd„  ■ cylinder  drag  coefficient  ■ drag  per  unit  length 

c 2aQ* 

®8 

cdc  • 1.2  for  laminar  crossfow;  cd(,  ■ .4)  for  turbulent  crossflow. 

Because  of  pressure  gradients  in  the  flow  field  in  which  the  store  is  immersed, 
buoyant  forces  will  also  appear. 


By  assuming  that  the  flow  in  planes  perpendicular  to  the  body  obeys  the  Laplace 
equation,  buoyant  force  terms  can  be  derived.  The  following  formulas  result. 


(CN)By 


dw[ 

dx] 


dx- 


(5.66) 


. 2n  f lt  2 dv®  J 

Cy  By  SRJ  a dxs  dxs 


dw! 


(Cm)By  * 4*ir£s  a2  as;  (x«.m ' V dx 


(5.67) 


(5.68) 


(Cn)By  * sfl ;jT  "3  a?  (xB.m  - xs)  dx 


5.2.2  Forces  on  store  lifting  surfaces 


(5.69) 


The  forces  on  lifting  surfaces  are  determined  by  assuming  that  the  angle  of 
attack  at  each  spanwlse  station  is  constant  and  given  by  the  value  at  the  quarter-chord 
point.  Then  reverse  flow  theorems  indicate  that  the  forces  and  moments  can  be  calcu- 
lated by  integrating,  across  the  span,  the  product  of  local  flow  angle  times  an  appropri- 
ate influence  function.  The  influence  functions  are  the  spanwlse  distributions  of  force 
or  moment  for  the  wing  in  a reversed  uniform  flow  field. 


First  it  is  necessary  to  obtain  the  component  of  velocity  distribution  normal 
to  the  lifting  surface  along  the  mean  quarter  chord.  The  velocities  in  the  store 
coordinate  system  ua,  vs,  wa  are  resolved  to  provide  the  components  normal  to  the  fins. 
The  lifting  surfaces  are  assumed  to  be  rotated  by  an  angle  with  respect  to  the  store 
coordinate  system,  as  shown  in  Fig.  5.15  for  a cruciform  arrangement. 

The  velocity  components  normal  to  surfaces  1 to  ft,  respectively,  are 


w1  * wa  cos  ♦f  + vs  sin  ♦f 
W2  ■ wa  cos  + va  sin 
vi  • -wa  sin  ♦f  + va  cos  ♦f 

v2  " -ws  sln  *f  + vs  cos  ♦f 


(5.70) 


Fins  Fi  and  Fj  are  the  rotated  "horizontal"  lifting  surfaces,  while  F3  and  F4  denote  the 
rotated  vertical  ones. 


The  calculations  are  simplified  by  defining  symmetric  and  unBymmetric  components 
of  the  spanwlse  angle  of  attack  distribution.  Thus  the  symmetric  component  is 


a 


B 


»i(yf)  + »»2  (-yf ) 


2V 


•s 


(5.71) 


and  the  unsymmetric  part  is 
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FIG.  5.15  COORDINATE  SYSTEMS  USED  IN  EMPENNAGE  FORCE  AND  MOMENT  CALCULATION 


Wi<yf)  - w2(-yf) 


(5.72) 


Then  the  total  force  along  the  z direction  (Including  body  carry-over)  in  the  fin  coordinate 
system  is 


[(cz)rh  " (cs 


>B]f ' * h “s(cc*>3  dyf 


where  from  slender  body  flow  for  the  fin-body  in  reverse  flow,  (ref.  60), 


(CCt)3  - 4 


/(Sh2yf2  - a^)(Sh2  - yf2) 


Sh2  yf2 


(5.74) 


(da  )h  1b  the  linear  theory  lift  curve  slope  of  the  horizontal  external  surfaces.  Joined 
together.  This  formula  uses  the  slender  body  theory  lift  distribution  but  corrects  the 
integrated  value  by  multiplying  by  the  ratio  of  linear  theory  lift  to  slender  body  lift. 

The  values  of  wi(y«.)  and  W2(-yf)  are  determined  from  the  previously  described 
flow  field  calculations  without  the  store  present.  The  store  Induces  a "Beskin"  upwash 
over  the  fin  panels  which  has  the  distribution 


wB(yf)  - w0  (5-75) 

w0  is  the  velocity , component  at  the  center  of  the  body  in  the  Zf  direction.  Thus, 
Instead  of  Eq.  (5.71)  for  yf  > a 

wi(yf)  + w,(-yf)  w 2 

' 2V«°  — ~ + £ ~2  (5-76> 

s “s  yf 

Similarly  the  force  along  the  yf  direction,  including  body  carry-over  and 
upwash,  is 


i 


! 
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where 


and 


(S&.J 


C(Cy)BV  - (C  )B]f  - I Bg(CCJl)j)  dzf 

<B„-a)2ya 


(CCt)j,  - n 


(sv2zf2  - a4)(sv3  - zf2) 


Bv2  *f2 


v2lZf)  + Vi(-zf)  v0a2 

8.  . -f— ± £-  + -2 ; Zf  > a 

* 2V-  ..  - 2 i - 


,dCL 


'»  V-8zf 


(5.77) 


(5.78) 


(5-79) 


Alao,  here  (sx^v  is  the  linear  theory  lift  slope  for  the  vertical  fins  and  v0  Is  the 


axso,  nere  itx- ;v  18  tne  linear  cneory  slope  lor  cne  ve: 

lateral  velocity  (along  the  yf  direction)  at  the  body  center 


The  forces  are  assumed  to  act  at  the  mean  quarter  chords  of  the  exposed  panels. 
Then  the  moment  coefficients  are 

[(Cm>BH  - <Cm)B]f  ■ - C(CZ)BH  - (Cz)B]f  (5.80) 

R 


C(Cn)Bv  “ " • C(Cy)Bv  “ Ir  (5»8l) 

where  lh  and  l are  the  distances  from  the  moment  center  to  the  fin  quarter  chord 
positions  (positive  for  the  fin  aft  of  the  moment  center);  and  lB  1b  the  reference 
dimension  for  pitching  and  yawing  moments. 


Rolling  moments  are  derived  for  the  planar  wing  case  and  the  cruciform  fin 
configuration  In  which  all  panels  are  of  identical  geometry.  The  procedure  Is  the  same 
as  for  the  other  forces  and  moments,  except  that  the  load  distribution  corresponds  to 
the  slender  body  theory  for  an  angle  of  attack  distribution  Increasing  linearly  along 
the  span.  The  results  are,  for  the  planar  case: 


In  this  formula, 


(ci)h 


<dS~)H 

(sh-' 


>H  f' 


‘au  (CC, 


>5  «*f 


(5.82) 


where  p Is  the  roll  rate, 
roll. 


wi(yf)  - w2(-yf)  pyf 


2V-. 


v._ 


(5.83) 


■s  -s 

and  the  last  term  in  Eq.  (5.83)  accounts  for  the  damping  In 


Slender  body  theory  for  a linearly  varying  angle  of  attack  (wing  with  linear 
twist)  gives 


(CC»)5  ■ (1  + f 008-1  —rj 


2ash 

— ^-r)(yf  ■ yr 

Sh2+a2  yf 


a£  t/(sh2-yf2)(3h2yf2-a4) 

* yr  1 


sh2yf2 


(yf2+a2)(Sh2-a2)^ 

r 


♦ jf  (yf  - — )2  cosh-1  [— ^5 — - — ^ — — ] 
8 yf  (yf2-a2)(Sh2+a2) 


For  cruciform  fins 


. (S£i)H  fSh 

(°Z)HV  “ " * r. ST  A (““  + ®u)(CCt)6  dyf 


(s-a)2l 

ou  Is  given  in  Eq.  (5.35),  and  Bu  is  obtained  from 

Vl(-Zf)  - v2(zf)  pzf 

Bu  - + — 

2V_  V. 


(5.84) 


(5.85) 


's 


's 


(5.86) 


»)V. 
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The  slender  body  result  for  the  cruciform,  derived  In  ref.  61  becomes 

(CC,)6  ■ [cos  29  tanh'1  (Bln  26  ) - cos  2y  tanh'1  (tan  29  )] 
1 0 » sin  2y  tan  2Y 


+ 4JL  [K(k:)  sin  *48 


- 2k1  cos  A1K(k1)Z(A1,k1)] 


(5.87) 


where 


r . 


Shz 


cos  26 


Sh2  (a1*  + yf  14 ) 

yf2  (a1*  + S^) 


(5.88) 


(5.89) 


cos  2y  ■ - — 

2R2 

ki  « sin  2y 

lUk^)  * complete  elliptic  Integral  of  the  first  kind 

K(kx)  •/ 

J ft 


IT  — 

7 dz 

/l-k^sln2  z 
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FCAj.kj)  » Incomplete  elliptic  Integral  of  the  first  kind 
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(5.90) 

(5-91) 

(5.92) 

(5.93) 

(5.9*4 

(5.95) 

(5.96) 

(5.97) 


At  yf  ■ a,  the  expression  becomes  Indeterminate  having  the  limit 

(CCfc)g  ■ [-  cos  2 y log(cos  2Y)+  .SKfkj)  sin  49]  (5-98) 

y*a 

5.3  Calculation  of  Store  Trajectory 

The  trajectory  of  the  store  Is  determined.  In  ref.  56,  by  Integrating  the  equa- 
tions of  motion.  The  trajectory  Is  determined  with  respect  to  the  airplane  (in  the  C, 
n,  c coordinate  system  fixed  to  the  fuselage).  The  parent  aircraft  Is  assumed  to  be 
flying  at  constant  velocity,  constant  angle  of  attack  and  constant  flight  path  angle 
with  respect  to  the  horizontal. 


5? 


C+[y  cos  (C,z)-z  cos  (;,y)]p+[z  oos  (C,x)-y  cos  (c,z)]q+[x  cob  (C,y)-y  cos  (C,x)]f 
Fx  . 

* [jr  + *(<i2  + p2)  - ypq  - zpr  - qz0  + ry0]  008  (c»*> 

+C^  - xpq  + y(p2  + r2)  - zqr  + pz0  - rx0]  cos  (c,y) 

Fz  , 

+C5T  - xpr  - yqr  + z(p2  + q2)  ♦ qx0  - py0]  cos  (t,z)  (5.101) 

The  rotational  equations  of  motion  are 

m[y  cos  (z,£)  - z cos  (y,?)]?  + m [y  cos  (z,n)  - z cos  (y,n)]ri 

+ m [y  cos  (z.O  - z cos  (y.O]?  + Ixxp  - Ixyq  - Ixzr 
« Mx  - rq(Izz  - Iyy)  + (q2  - r2)Iyz  + p(qlxz  -rlxy) 

- m [7(py0  - qx0)  - z(rx0  - pz0)]  (5.102) 

m [z  cos  (x,t)  - x cos  (z,C)]£  + m [z  cos  (x,n)  - x cos  (z,n)]n 


+ m [z  cos  (x,c)  - x cos  (z,t)]c  - IxyP  + Iyyq  - Iyzr 
-My  - rp(Ixx  - Izz)  + (r2  - p2)Ixz  + q(rlxy  - plyz) 

- m [z(qz0  - ry0)  - x(py0  - qx0)]  (5.103) 
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m [x  cos  (y,$)  - y cos  (x,e)](  + m [x  cos  (y,n)  - y cos  (x,n)]H 
+ m lx  cos  (y ,; ) - y cos  (x,;)];  - Ixzp  - Iyz<i  + Izzr 
■ Mz  - pqUyy  - Ixx)  + (P2  - q2)IXy  + r(plyz  - qlxz) 

- ra  [x(rx0  - pz0)  - y (qzD  - ry0)]  (5 . 104 ) 

The  quantities  (,  n,  (,  p,  q,  r are  determined  by  Integrating  the  six  equations 
of  motion,  €,  Hi  (are  the  coordinates  of  the  moment  center  of  the  store  with  respect 
to  the  parent  aircraft,  and  p,  q,  r are  the  angular  velocities  of  the  store  about  the 
store  axes  (x,  y,  z)  illustrated  in  Pig.  5.13. 


The 

which  locate 
coordinates. 

equations  account  for  mass  asymmetry  through  the  coordinates  x,  y,  z 
the  center  of  mass  with  respect  to  the  moment  center  in  the  store  body 

The 

direction  cosines 

relate  the  store  axes  to  the  parent  aircraft  axes. 
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(5.105) 

As  indicated  in  Pig.  5. 11*,  the  store  axes  are  related  by  rotations  of  yaw,  ill 
pitch,  6,  and  roll,  $ (in  that  order)  to  the  aircraft  axes.  These  angles  are  obtained 
by  time  integration  of  the  stores  angular  rates. 

^ „ (q  sin  ♦ + r cos  ») 
cose 

6 ■ q cos  ♦ - r sin  ♦ 

$ ■ p + q sin  $ tan e + r cos  e tan  9 (5.106) 

The  integrations  are  performed  numerically  starting  with  some  prescribed  initial  position. 

The  quantities  x0,  y0,  and  zQ  are  the  velocity  component-  of  the  store  moment 
center  in  the  store-fixed  coordinates.  These  are  related  through  the  acceleration 
equations 

3c0  - c cos  (x,()  + n OOS  (x,n)  + ii  COS  (x,() 

?0  - C cos  (y,£)  + ri  cos  (y,n)  + t cos  (y,c) 

z0  • 5 cos  (z,£)  + ii  cos  (z,n)  + C cos  (z,C)  (5.107) 

The  forces  and  moments  are  obtained  by  adding  all  of  the  aerodynamic  contribu- 
tions to  the  gravitational  forceB. 
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| *X  * <^X  - 1-5SHCA 

I 

r 

Py  - mgy  + q»gSRCy 

Pz  ’ (5.108) 


Mx  - m(gzy  - gyz)  ♦ q.|JSRtRClt 
My  - m(gx  z - Bzx)  + q-BSRtRCm 

Mz  « m(gyx  - gxy)  + 9.8sRiRcn  (5.109) 

Hera  m Is  the  mass  of  the  store,  la  a prescribed  axial  force  coefficient,  and  the 
other  aerodynamic  coefficients  are  the  total  contributions  to  the  forces  and  moments 
of  all  components  of  the  store  configuration. 

The  components  of  gravity  are  transferred  fr.jm  the  parent  aircraft  coordinates 
to  the  store  coordinates  to  give 


gx  » -g  cos  e cos  4 sin  (of+yf)  - 45  sin  8 cos  (af+Yf) 


I 

i 


i 


1 


gy  ■ -g(sln  4 sin  8 cos  4 - cos  4 sin  4)  sin  (of+Yf)  + g sin  4 cos  8 cos  (af+Yf) 

i 

gz  ■ -g(cos  4 sin  8 cos  4 + sin  4 sin  4)  + g cos  4 cos  8 cos  (af+Yf)  (5.110) 


ar  Is  the  angle  of  attack  of  the  parent  aircraft  and  yf  is  the  angle  between  the  air- 
plane's flight  path  and  the  horizontal  (positive  when  the  airplane  Is  climbing). 

Finally,  the  moments  and  products  of  Inertia  are  defined  by 


Ixx  ’ /(y2+z2)dm 


Iyy  ■ /(x2+z2)dm 

Izz  “ /(x2+yz)dm 


IXy  - [x y dm 


Ixz  ■ /xz  dm 


Iyz  ■ /yz  dm  (5.111) 


Computer  programs  for  carrying  out  all  of  these  computations  are  described  in 
Vol.  II  of  ref.  58. 

5.*i  Comparisons  With  Experiment 

5.0.1  Plow  field 

Comparisons  of  the  computational  procedures  with  experiments  are  reported  in  j 

ref.  56.  A few  illustrations  from  that  report  will  be  shown  here.  Examples  of  flow 
field,  forces  and  trajectories  will  give  an  Indication  of  the  accuracy  of  the  method. 

) 

A wing  fuselage  model  is  shown  in  Pig.  5.16.  A store  location  on  this  config- 
uration Is  shown  in  Fig.  5.17,  and  store  details  In  Fig.  5*18.  In  Fig.  5-19  sldewash 
and  upwash  distributions  along  the  store  centerline  (without  the  store  present)  are 
compared  with  the  theory  for  the  aircraft  at  M.  ■ .25  and  6-degree  angle  of  attack. 

The  agreement  Is  quite  good.  The  analysis  accurately  predicts  the  change  in  upwash 
due  to  a pylon,  but  is  not  quite  so  effective  on  the  sldewash  increment. 


P 


PYLON 


PIG.  5.20  EFFECT  OF  PYLON  ON  LOAD  DISTRIBUTION 
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5.^.2  Force  Distribution 

The  force  distributions  on  the  store  are  shown  In  Figs,  5. 20. a and  5.20.b. 

Again  the  theory  Is  quite  accurate. 

5.4.3  Tra.lectorv 

The  measured  Initial  trajectory  of  the  store  (by  the  captive-store  teat 
technique),  is  compared  with  calculations  In  Fig.  5.21.  In  this  case  the  Mach  number 
Is  increased  to  .4  while  the  angle  of  attack  Is  4 degrees.  The  pylon  is  included.  The 
displacements  of  the  midpoint  of  the  store  are  quite  accurately  predicted  by  the 
calculation.  The  pitch  angle  is  somewhat  off,  and  would  presumably  result  in  some 
divergence  of  the  vertical  displacement  at  later  times.  The  gravitational  force  is  a 
significant  fraction  of  the  total  however,  so  that  the  trajectory  is  not  too  sensitive 
to  store  angle.  The  store  without  fins  would  be  unstable.  A stable  store  would  probably 
not  reach  such  high  angles  of  attack  under  these  release  conditions. 

The  high  predicted  value  of  pitch  angle  is  consistent  with  the  force  distribu- 
tion shown  in  Fig.  5. 20. a.  The  theory  shows  a little  too  much  normal  force  on  the  nose 
and  a high  negative  force  on  the  rear.  The  extra  pitching  moment  implied  by  these 
discrepancies  would  result  in  higher  pitch  angle  predictions  than  would  be  measured  in 
the  captive  trajectory  experiment.  However  the  magnitude  of  the  error  is  larger  than 
would  be  expected. 

On  the  whole,  the  trajectory  predictions  are  very  good  for  evaluation  of 
release  safety  and  could  give  good  Impact  point  predictions  for  stable  stores. 


CAPTIVE  STORE 
NO  DAMPING 


THEORY 
NO  DAMPING 


FIG.  5.21  COMPARISON  BETWEEN  CALCULATED  TRAJECTORY  AND  CAPTIVE-STORE  TRAJECTORY 
OF  STORE  S0  RELEASED  AT  ONE-THIRD  SEMI-SPAN  LOCATION;  M.  = 0.4 
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5.5  Supersonic  Plow  Field 

5.5.1  Introduction 

The  calculation  of  the  trajectory  of  a store  released  at  supersonic  speed 
proceeds  in  exactly  the  same  manner  as  the  subsonic  case.  The  flow  field  distribution 
at  the  store  location  is  first  determined;  from  this  distribution  the  forces  and  momentB 
acting  on  the  store  are  computed,  and  then  the  trajectory  is  obtained  by  integrating  the 
resulting  accelerations. 

Unfortunately  the  flow-field  calculation  is  complicated  by  the  presence  of 
shock  waves  in  the  flow.  Hence,  while  linear  methods  are  available  for  predicting 
supersonic  flow  fields,  they  are  not  sufficiently  accurate  in  locating  shock  waves  nor 
in  predicting  the  increments  in  the  flow-field  parameters  across  the  shocks.  Consequently 
some  augmentation  of  linear  methods  is  required. 

The  procedure  described  here  follows  the  method  of  F.  D.  Fernandes;  it  is 
described  in  more  detail  in  refs.  62  and  63.  The  subsonic  analysis,  ref.  62,  is  very 
similar  to  the  method  of  Goodwin,  Dillenius,  and  Nielsen  described  in  Sections  5.1  and 
5.2  (refs.  56-58).  The  computer  programs  developed  by  Fernandes  are  described  in 
ref.  66. 


One  difference  between  the  method  of  Fernandes  and  that  of  Goodwin,  et  al,  is 
in  the  method  of  computing  the  forces  on  the  store.  Fernandes  determines  Increments 
due  only  to  the  Interference  flow  field  and  adds  them  to  the  store-alone  forces  in  a 
uniform  flow.  Goodwin,  et  al,  calculate  the  forces  based  on  the  entire  local  flow 
field.  There  are  small  differences  in  the  resulting  aerodynamic  coefficients  when 
vortex  separation,  for  example,  introduces  nonlinearities. 

After  shock  waves  are  located  by  means  of  a geometric  transformation,  the  flow 
field  about  aircraft  components  is  determined  by  linear  theory.  Then  the  load  distribu- 
tion on  the  store  is  calculated  by  dividing  the  store  into  axial  sections,  and  determining 
the  load  on  the  sections  by  linear  crossflow  methods.  The  procedure  calculates  loading 
on  the  body  due  to  local  crossflow,  axial  rate  of  change  of  crossflow,  and  buoyancy. 

Fin  loads  due  to  crossflow,  and  chordwise  and  spanwise  rates  of  change  of  crossflow  are 
taken  into  account.  Loads  due  to  vortex  separation  and  fin-body  interference  are 
included. 


5.5.2  Location  of  shock  waves 
5. 5. 2.1  Attached  shock  waves 

To  account  for  the  finite  strength  and  curvature  of  bow  or  leading  edge  shock 
waves,  the  actual  body  is  transformed  into  an  "equivalent  body"  as  indicated  in  Fig.  5.22. 
Here,  at  some  lateral  distance  from  the  body,  a Mach  line  1b  drawn  through  the  actual 
shock  wave  to  the  body  axis.  Then  an  enlarged  body  is  constructed  by  geometrically 
scaling  up  the  actual  body  so  that  the  noBe  of  the  transformed  body  coincides  with  the 
origin  of  the  Mach  line,  while  the  shoulder  station  remains  fixed.  As  in  the  figure, 
the  local  radius  of  the  scaled-up  body  is  related  to  that  of  the  actual  body  by  the 
formula 

R'  . r*  MlM  (5.112) 

X« 


FIG.  5.22  EQUIVALENT  BODY  CONCEPT  FOR  IMPROVEMENT  OF  LINEAR  THEORY 
FLOW  FIELD  CALCULATIONS  IN  SUPERSONIC  FLOW 
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where  XS  * distance  from  body  shoulder  to  nose,  and  XD  » forward  shift  required  to  make 
the  nose  coincide  with  the  Mach  line. 

Then  the  flow  field  at  the  chosen  radial  distance  is  found  from  linear  theory 
ror  the  transformed  body.  The  scale  of  the  transformation  will  vary  with  distance  from 
the  body  since  the  Mach  line  will  be  traced  back  from  different  points  on  the  bow  shock 
wave  of  the  real  body. 

The  problem  Is  to  find  the  shock  shape  and  location.  For  this  purpose  the 
procedures  of  Love  and  Long  are  employed.  The  shape  of  attached  shock  waves  Is  given  In 
rer.  65,  while  detached  shapes  are  treated  in  ref.  66.  From  ref.  65,  the  shape  of  the 
attached  shock  wave  is  given  by  the  formula 


yA  - tan  e tn(l  + x*)  + tan  pCx*  - ln(l  + Xi)] 


(5.113) 


xi  and  yj  are  the  coordinates  of  points  on  the  shock  wave,  c Is  the  shock  wave  angle 
corresponding  to  the  Initial  slope  of  the  body  nose,  6,  and  the  Mach  number  of  the  flow, 
M.  For  the  two-dimensional  case,  such  as  a wing  or  pylon  section,  e is  determined  from 
oblique  shock  theory;  while  for  bodies  of  revolution  the  solution  for  the  tangent  cone 
is  used . The  functions  may  be  found  in  ref,  67,  u is  the  Mach  angle  given  by 


sln"1  I 


(5. lH) 


Plots  ol  shock  v,'&ve  angle  £ as  functions  of  Mach  number  and  turning  angle  are 
shown  in  Figs.  5-23  snd  5*24  for  the  two-  and  three-dimensional  cases , respectively. 

Although  Eq.  (5.113)  gives  the  shape  of  the  attached  shock  wave,  the  scale  must 
be  tied  to  the  .iize  of  the  nose  of  the  body.  This  relation  is  given  by 


({-) 


KX< 


(f-)  - Kyi 


(5.115) 


!c  is  the  length  of  the  nose  of  a circular  arc  fitted  to  the  forward  portion  of  the  actual 
nose  shape.  K Is  a scale  factor  determined  by  fitting  Eq.  (5-113)  to  shock  shapes  for 
circular  arc  noses.  K Is  a function  of  Mach  number  and  turning  angle.  The  results  of 
the  curve  fits  are  tabulated  below  and  plotted  in  Figs.  5.25  and  5.26. 

TABLE  5.1  VALUES  OF  K 
(a)  for  two-dimensional  nose  shapes 

Values  of  K for  5 ■ 


m 

5° 

O 

O 

15° 

20° 

25° 

306 

350 

1.5 

8.0 

5.0 

_ 

_ 

_ 

m , 

2.0 

8.0 

5.0 

3.0 

1.5 

- 

- 

- 

2.6 

8.0 

5.0 

3.0 

2.0 

1.15 

0.50 

- 

3.3 

8.0 

5.0 

3.0 

2.0 

1.50 

1.10 

0.H5 

<4.0 

8.0 

5.0 

3.0 

2.0 

1.50 

1.10 

.70 

5.0 

8.0 

5.0 

3.0 

2.0 

1.55 

1.20 

.85 

H0o 


0.A3 


(b)  For  axlsymmetrlc  nose  shapes 
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/ 

FIG.  5.27  COMPARISON  OF  PREDICTED  AND  EXPERIMENTAL  SHOCK  FOR  AXISYMMETRIC 
NOSE  WITH  CONTOUR  OF  THE  FORM  y = a - bx2  + cx3  - dx1*  (WITH  ORIGIN 
FOR  NOSE  CONTOUR  AT  BASE  OF  NOSE).  M = 1.62;  6 = 28° 


A comparison  of  theory  and 
experiment  for  a particular  axlsymmetrlc 
nose  Is  shown  In  Fig.  5-27. 


5. 5. 2. 2 Detachment  distance 

For  detached  shocks,  the 
detachment  distance  is  required  aB  well 
as  the  shape  of  the  shock  wave.  The 
formula  for  detachment  distance 
employed  In  ref.  66  is 

if  ■ .5C  cot  «det  (5.116) 

x'  Is  the  detachment  distance  and  d' 
the  diameter  or  thickness  of  the  body 
nose  as  shown  In  Fig.  5.28.  Sdey  is 
the  cone  or  wedge  angle  that  would 
Just  cause  shock  detachment  at  the 
stream  Mach  number. 


Fernandes  (ref.  63)  found  that 
this  procedure  gives  poor  results  for 

some  ogive  cylinder  bodies.  Therefore  FIG.  5.28  DETACHED  SHOCK  GEOMETRY 

he  replaces  ogive  cylinders  by  "equiva- 
lent" cone  cylinders.  The  modification 
of  the  body  Is  shown  in  Fig.  5.29. 


FIG.  5.29  MODIFIED  GEOMETRY  FOR  OGIVE-CYLINDER  BODIES 
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A cone  is  constructed  through  the  apex  of  the  nose  and  the  point  at  which  the 
slope  of  the  nose  is  equal  to  “det-  The  cone-cylinder  thus  formed  is  analyzed  by  the 
procedure  for  cone-cylinder  bodies  to  be  described  a little  further  on. 

(Flat-faced  bodies) 

For  a cone  or  wedge  with  nose  angle  equal  to  {det,  the  constant  C in  Eq.  (5.116) 
would  be  1.  For  other  shapes,  C has  slightly  smaller  values.  For  axisymmetrlc  bodies 
with  flat  faces,  by  comparison  with  experiment, 


For  two-dimensional  flat-faced  bodies, 


(5.117) 


(5.118) 


Figure  5.30  (taken  from  ref.  66)  shows  a comparison  of  Eq.  (5.116)  with  experimental 
data  from  refs.  68-78  for  a range  of  supersonic  Mach  numbers. 

(Circular-faced  bodies) 

For  bodies  with  complete  semicircular  noses,  such  as  spheres,  hemisphere- 
cylinders,  or  the  two-dimensional  counterparts  (circular  cylinders,  or  semicircular- 
nosed  flat  plates),  the  appropriate  values  of  the  coefficient,  C«Cc , vary  with  Mach 
number.  This  functional  dependence  is  shown  in  Fig.  5.31  for  two-dimensional  and 
axisymmetrlc  bodies. 

(Cone-cylinders  and  wedge  slabs) 

For  cone-cylinder  bodies,  or  the  analogous  wedge-slab  two-dlmenalonal  configura- 
tion, the  data  for  detachment  distance  is  fitted  by  the  formula 


' -5  cot  «det  + y 


(5.119) 


where  y is  the  negative  solution  of  the  quadratic  expression 


<£i»!  - * c15#  * h »-*>2  * n ■ » 


(5.120 


2N  - Iq 


(5.321) 


N ■ .5  cot  5det  (C90-l) 


(5.122) 


b2  - ( b 4 N ) 2 


(5.123) 


1 ■ - - (radians) 


(5 . 121) ) 


- .5  csc^  S 


(5.125) 


x * 60  - 5det  (radians) 


(5.126) 


and  S0  is  the  cone  or  wedge  semi-apex  angle. 
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SEMIAPEX  ANGLE,  4„,  DEG 

FIG.  5.J2  VARIATION  IN  DETACHMENT  DISTANCE  WITH  SEMI-APEX  ANGLE  FOR  CONE- 
CYLINDERS  AND  WEDGE-SLABS  AT  CONSTANT  MACH  NUMBER 


Comparisons  of  this  formula  with  data  from  ref.  79  are  shown  In  Fig.  5.32. 

(Cut  spheres  and  cylinders) 

From  ref.  66,  the  variation  of  detachment  distance  for  cut  spheres  (or  cut 
cylinders  In  the  two-dimensional  case)  Is  given  in  the  form 


xj. 

d' 


•5C  cot  «det 


(5.127) 


A cut  cylinder  or  sphere  Is  Illustrated  In  Fig.  5*33.  In  this  equation,  C Is  obtained 
from 


C = C0 


2 


sdet 


< 0 < 


? 

(5.128) 


C + y'  0 < 6 < 


- - ? 


- 6 


det 


where  Cc  Is  the  value  for 
complete  spheres  or  cylinders 
obtained  from  Fig.  5-31. 

As  in  the  cone- 
cylinder  case,  an  elliptic 
variation  of  C with  60  (where 

60  * \ -0)  produces  the 

equation 


SHOCK  WAVE 


FIG.  5.33  CUT  CYLINDER  OR  SPHERE  WITH 
DETACHED  SHOCK  WAVE 


r2-  - 1]  - 0 


(5.129) 


q (b+N ) 


(5.130) 


211  - l q 


(5.131) 


.5  cot  6det 


(5.132) 


(radians ) 


11  ■ Cgg  - Cc 


(5.131!) 


p ~ 0 - 6det 


(radians) 


(5.135) 


The  function  needed  to  define  q In  Eq.  (5.132),  Is  obtained  by  matching  the 
slope  at  C * Cc  of  the  ellipse  through  the  point  C ■ 1 at  S0  » Sdet  to  the  slope  at  C » C0 
of  the  ellipse  defined  by  Eq.  (5.129).  The  equation  of  the  ellipse  through  C » 1 at 
40  - 6det  and  having  Its  minimum  at  60  « 90  degrees,  Is  given  by 


1 y2  . C2_(b+N)  ]y  +[Cb+tjll  + [x-lf  . x]  . 0 


(5.136) 


a,  b,  1,  x are  defined  by  Eqs.  (5.130),  (5.13D,  (5.133)  and  (5.135),  respectively, 
but  now 


q * - -5  esc  6 


(5.137) 


fJ  - .5  cot  6det  (C90  - 1) 


(5.138) 


At  the  desired  point,  y * .5  cot  S*-*  (C„  - 1)  and  x ■ the  corresnondlng  value  from 
(5.136) 


U . . b2  r x-L 
^ Ly  - (b+N) 


(5.139) 


This  slope  defines  q in  Eq.  (5.132)  so  that  y'  can  be  found  as  a function  of  x (or  6) 
from  Eq.  (5.129)  and  C from  Eq.  (5.128) 
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A comparison  of  predicted  and  measured  shock  detachment  distances  for  cut 
spheres  (attached  to  cylindrical  afterbodies)  is  shown  in  Fig,  5 . 34 . The  plot  shows 
detachment  distance  measured  by  the  parameter  shown  on  the  figure,  which  can  be 

related  to  the  predicted  values  of  54-  through  the  relations 

d ' 


2 tan  6 2 sin  6 


0 < 8 < 1 - 6 

“ “2 


det 


b . Bln  <?  - w 

D sin  6 


d' 


tan  - 4det> 


2 sin  6 


\ “ 4det  < 6 < 90° 

(5.140) 


5. 5. 2. 3 Shape  of  detached  shocks 


The  shape  of  the  detached  shock  wave  is  determined  in  ref.  66  by  a modification 
of  Moeckel's  method  (ref.  80).  The  geometrical  arrangement  is  shown  In  Fig.  5.35. 

The  angle  0,  defining  the  control  line  is  determined  as  a function  of  Mach  number  by 
correlating  the  results  for  known  shock  shapes.  The  functions  for  spheres  and  cylinders 
are  shown  in  Fig.  5-36  along  with  functions  used  in  other  analyses. 


FIG.  5.34  PREDICTION  OF  RESULTS  AT  M«  = 3.55  FOR  A SPHERE  THAT  IS 
EFFECTIVELY  CUT  (DIAMETER  OF  ACTUAL  MODELS  HELD  CONSTANT 
AND  RADIUS  OF  NOSE  VARIED) 
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FIG.  5.35 


GENERAL  FEATURES  OF  THE  METHOD  OF  MOECKEL 
ADOPTED  IN  PRESENT  ANALYSIS 


FREE-STREAM  MACH  NUMBER,  Me, 
(«)  SPHERE 


AT  M„, 


{MAXIMUM  7,3 
FROM  C(pv 
TANGENTIAL"^*] 


VARIATION  OBTAINED 
IN  PRESENT  ANALYSIS 
(FROM  AVERAGE 
VALUES  OF  r,) 


1 2 3 4 5 6 

FREE-STREAM  MACH  NUMBER,  Moo 

(b)  CIRCULAR  CYLINDER 
(TWO-DIMENSIONAL) 


FIG.  5.36  COMPARISON  OF  VARIATION  OF  r,  WITH  MACH  NUMBER  OBTAINED 

IN  PRESENT  ANALYSIS  WITH  THOSE  OBTAINED  IN  OTHER  ANALYSES 
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With  this  point  established,  the  shock  shape  is  given  by 


i_  . I V(— )2  - (^)2 

a • 8 ’ld';  ld,; 


(5.1*11) 


The  distances  x and  y are  measured  from  the  origin  of  the  asymptote  ae  shown  in  Fig.  5.35* 

8«  /fo2-l;  d'  is  the  diameter  or  width  of  the  body  at  the  point  of  tangency  of  the  cone 
or  wedge  of  incipient  detachment. 


The  vertex  of  the  shock  is  given  by  the  formula  (from  Moeckel), 

x' 

(r> 


Xo 

d' 


8 / B2  tan2  e„  - 1 & + 


2 2 

8 tan^  t. 


8 


/ B2  tan2  c „ 


(5.1*t2) 


- 1 + tan  n 


where  e is  the  shock  wave  angle  that  will  Just  give  sonic  flow  behind  the  shock,  x'  is 
the  shock  detachment  distance  determined  previously  for  noses  of  various  shapes. 


Some  comparisons  with  experiment  are  shown  in  Fig.  5.37  taken  from  ref.  66. 
In  these  plots,  D is  the  diameter  of  the  cylindrical  body  while  F is  the  horizontal 
distance  measured  from  the  front  of  the  flat-faced  cylinder  or  the  center  of  the 
hemispherical-nosed  body. 

5.5*3  Supersonic  flow  field  determination 

5. 5* 3.1  Wing  and  pylon  representation 

With  the  Bhock  shape  established  as  indicated  in  the  previous  section,  the 
supersonic  flow  field  can  now  be  obtained  by  application  of  linear  theory  to  a body 
transformed  to  make  its  shock  wave  match  the  more  accurately  determined  location. 


FIG.  5.37  PREDICTION  OF  SHOCK  SHAPE  AND  LOCATION  FOR 
AXISYMMETRIC  NOSES  AT  M„  = 3.55 
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The  linear  employed  by  Fernandes  for  a wing  or  pylon  in  a crossflow  uses  the 
"Mach  box"  method  of  Zartarlan  (ref.  81).  The  perturbation  velocity  potential  due  to 
a source  distribution  is  given  by 


$ (x,j 


Mach 

forecone 


w(x1,y1)dx1dy1 
/(x-xj)2  - 62(y-y1) 


2 


(5.143) 


w(x1,y,)  is  the  distribution  of  downwash  velocity  in  the  region  of  the  Mach  forecone 
ahead  of  the  point  x,  y. 


In  the  "Mach  box"  method,  w(xT,yn)  is  assumed  to  be  constant  in  rectangular 
boxes  whose  diagonals  are  Mach  lines.  With  w constant,  the  integration  can  be  carried 
out  and  the  velocity  potential  at  x,y  can  be  found  for  a Mach  box  centered  at  any  posi- 
tion in  the  foi'econe.  By  shaping  the  rectangles  according  to  Mach  lines,  the  velocity 
potential  function  becomes  independent  of  Mach  number.  The  wing  planform  and  the 
adjacent  disturbed  region  of  the  flow  field  are  subdivided  into  Mach  boxes  as  shown  in 
Fig.  5.38  for  two  sample  configurations. 


The  numbers  in  the  boxes  designate  the  order  in  which  they  are  treated.  Since 
w is  prescribed  in  box  number  1 by  the  downwash  distributicn,  and  with  no  other  boxes 
in  its  foreoone , $ can  be  obtained  explicitly.  In  boxes  on  the  wing,  the  downwash  is 
known.  In  boxes  in  the  "diaphragm"  region,  the  pressure  coefficient  is  zero  since  the 
diaphragm  is  a region  of  unknown  downwash  which  deflects  to  balance  the  preB3ure. 


The  pressure  coefficient,  by  linear  theory,  is  proportional  to  the  perturba- 
tion velocity  in  the  streamwise  direction 


C - _ lu  , _ 2_  li 
U„  U»  3x 


(5.144) 


By  differentiation  of  the  velocity  potential  for  a Mach  box  and  summation  over  all  Mach 
boxes,  in  the  forecone  of  a given  point,  the  perturbation  velocities  can  be  determined. 


U1 


E 

in  l^ach 
forecone 


(5.145) 


R(1,J)  is  called  the  aerodynamic  Influence  coefficient,  and  is  a function  of  the  location 
of  Mach  box  J with  respect  to  the  point  at  i.  R(1,J)  is  the  derivative  of  the  velocity 
potential  due  to  the  Jth  Mach  box  for  unit  downwash. 


R(1,J) 


1 y - yi 
x |y  _ yxl 


[}yi  ' yjl]s 

^(Xi-Xj) 


(5.146) 


MACH  LINES 


FIG.  5.38  ORDER  OF  COMPUTING  MACH  BOX  PROPERTIES  (U  OR  W)  ON  WING  AND  DIAPHRAGM 
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where  the  subscript  f implies  that  the  function  is  evaluated  over  the  area  of  the  Mach 
box.  That  is,  the  sides  of  the  Mach  box  form  the  limits  of  Integration,  while  the 
function  is  the  integral  of  Eq.  (5.193)  for  w(xi,yi)  ■ 1. 

For  a given  subdivision  of  the  Mach  forecone  into  boxes,  R(1,J)  is  a function 
only  of  the  position  of  the  box.  Values  are  shown  in  Fig.  5.39. 

From  the  condition  that  Cp  ■ 0 on  the  diaphragm,  the  values  of  w (downwash) 
on  the  diaphragm  boxes  can  be  obtalhed  by  proceeding  along  the  boxes  in  order. 

5. 5. 3. 2 Flow  field  due  to  wing  and  pylon 

Once  the  values  of  w have  been  assigned  to  all  of  the  Mach  boxes  on  the  wing 
and  diaphragm,  the  velocity  potential  and  hence  the  perturbation  velocities  can  be 
determined  at  any  point  in  the  flow  field.  The  region  of  influence  is  illustrated  in 
Fig.  5-90. 


The  velocity  potential  at  the  point  x,  y,  z,  due  to  a Mach  box  in  the  z * zi 

plane  is 


f.ff 


dxjdyj 


/7*. 


•':i)"  - (z-z-^2  - (y-yj)2 


(5.197) 


where  is  the  area  of  the  Mach  box  (or  any  rectangular  area  over  which  w is  assumed 
to  be  constant). 
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FIG.  5.39  VALUES  OF  AERODYNAMIC  INFLUENCE  COEFFICIENTS  CAIC’S) 
VERSUS  BOX  LOCATION  IN  MACH  FORECONE 


MACH  BOX  IN 
THE  z.  PLANE 


FIG.  5.90  REGION  OF  X-Y  PLANE  INFLUENCING  THE  POINT  X,Y,Z 
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Carrying  out  the  Integrations  and  differentiating  to  obtain  the  velocity 
components  leads  to  the  results 

Case  I 


u(x,y,z)  - 0 


ttB 


v(x,y,z)  • Qv 

ti  e * 


wUi.y-i.zi) 

w(x,y,z)  - ^ Qz 


(5.188) 


The  functions  Gx,  Gy,  and  Gz  arise  from  the  Integration  and  differentiation  of  Eq.  (5-187), 
and  are  given  belovr 
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(5.189) 
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If  the  rectangle  lies  on  the  border  of  the  Mach  forecone,  then  some  of  the 
functions  will  not  be  real.  In  that  case  Xp  Is  taken  to  bs  the  most  forward  part  of  the 
forecone  that  lies  Inside  the  rectangle  (see  Fig.  5.81).  Then  the  derivatives  give 
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FIG.  5.81  POSITION  OF  MACH  BOXES  WITH  RESPECT 
TO  FORECONE 
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v(x,y,z) 
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(5.150) 


The  velocity  field  due  to  a distribution  of  rectangular  source  elements  Is 
obtained  by  summing  the  contributions  from  all  of  the  sources  within  the  Mach  forecone 
from  the  point  x,y,z. 


5.5. 3. 3 Flow  field  due  to  fuselage  nose 


The  flow  field  about  the  fuselage,  or  other  body,  is  determined  by  a Karman- 
Moore  procedure.  A description  of  the  theory  may  be  found  in  ref.  82.  A brief  summary 
will  be  given  here. 


The  fuselage  is  represented  as  a pointed  body  of  revolution.  For  the  body  at 
zero  angle  of  attack,  a line  distribution  of  sources  Is  placed  along  the  axis.  Assume 
that  the  sources  are  adjusted  at  the  points  x^,  x2,  x^.  • -as  shown  In  Fig.  5.42 


FIG.  5.42  SOURCE  POINTS  FOR  CIRCULAR  FUSELAGE 


The  perturbation  velocities  Induced  by  a line  of  sources  of  strength  gx  are,  from  linear 
theory  for  axlsymmetrlc  flow. 
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and 


(5.152) 


where  the  source  strength  varies  as  gx.  The  boundary  condition  requires 


vr  . dR 
U»+u  dx 


(5.153) 


The  first  source  distribution,  starting  at  x^,  produces  a conical-tipped  nose. 
The  slope  at  will  be 
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(5.154) 


Hence,  to  match  the  known  slope  of  the  body  at  Pj 


El 


dR 

3x 


air 

37 


ll  cosh-1  (i%)  + 9/(J|_,2  - 


(5.155) 


With  this  value  established,  the  change  in  source  strength  distribution  at  x2  can  be 
determined  to  make  the  slope  of  the  flow  at  P2  match  that  of  the  body  there.  The  formula 
for  the  slope  g^  is 
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(5.156) 
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To  account  for  angle  of  attack,  the  line  sources  are  supplemented  by  line 
doublets.  The  flow  field  due  to  a line  doublet  whose  strength  Increases  linearly  with 
slope  c results  In  velocity  components  (ref.  83,  for  example). 
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(5.159) 


where  e is  the  angle  around  the  body  measured  from  the  windward  element. 
The  boundary  condition  on  the  body  surface  requires  that 
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(5.160) 


Again  changes  In  doublet  strength  slope  can  be  adjusted  to  match  the  boundary  condition 
at  selected  axial  stations.  The  resulting  formula  for  the  slope  c^  Is 


, dR  | 
1 + 3x| 
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(5.161) 


The  flow  field  around  the  body  Is  determined  by  the  velocities  obtained  from 
Eqs.  (5.151),  (5.152),  (5.157),  (5.158),  and  (5.159).  The  source  strength  distributions 
are  determined  by  Eqs.  (5.156)  and  (5.161).  The  pressure  on  the  body  can  be  found  from 
the  linear  theory  formula 
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£u  ,'rii 

u 


(5.162) 


A comparison  of  measured  and  calculated  pressure  distributions  on  a body  at 
zero  angle  of  attack  is  shown  In  Fig.  5.43  (from  ref.  63).  The  good  agreement  13  an 


FIG.  5.43  COMPARISON  OF  THEORY  WITH  EXPERIMENT  FOR  PRESSURE 
COEFFICIENT  ON  AXISYMMETRIC  NOSE 
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Indication  that  the  flow  field  calculation  is  also  accurate  close  to  the  body  surface. 
Far  from  the  body  axis,  the  linear  theory  Is  not  so  accurate,  hence  the  shock  wave 
correction  Is  applied  to  transform  the  body,  as  described  earlier.  The  body  flow  field 
calculation  cannot  be  applied  when  the  first  Mach  line  passes  Inside  the  body  contour. 
Thus  blunter  bodies  can  only  be  treated  by  making  additional  approximations. 

5.6  Forces  on  Store  In  honunlform  Supersonic  Flow 

5.6.1  Store  body 

The  force  distribution  on  a store  in  a nonuniform  flow  field  is  calculated 
In  ref.  63  by  subdividing  the  store  into  axial  sections  having  nearly  uniform  angle  of 
attack.  The  normal  force  per  unit  length  on  each  section  Is  assumed  to  have  the  form 


■'N, 


“ 33T  + CN  57 


S(x) 


For  example,  by  slender  body  theory 

CN(x)  - 2a  -p- 

where  S Is  the  cro.  s-sectlonal  area  at  station  x and  „R2  Is  the  base  area. 
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Then 
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Hence  In  this  case 
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In  general  the  functions  and  CN  can  be  obtained  from  experiment  or  linear 
theory  as  used  In  the  supersonic  flow  field  determination.  The  values  chosen  for  Cfj  and 
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and  are  those  at  the  middle  of  each  axial  section.  The  average  value  of  the  angle 

the  angles  at  the  c 

6ui  + <*1-1  + 31 +1 


of  attack  Is  obtained  by  connecting  the  angles  at  the  centers  of  the  axial  sections  by 
straight-line  variations.  Then 


al 


(5.166) 


The  slope  ^ is  also  obtained  from  the  average  angles  found  by  Eq.  (5.166). 

In  Fernandes'  program,  nonlinear  effects  due  to  body  vortex  separation  are 
accounted  for  only  by  linearizing  about  the  lift  curve  slope  at  some  finite  angle  of 
attack  instead  of  the  zero  lift  slope. 

Buoyant  forces  on  the  store  are  determined  from  pressure  differences  in  the 
flow  field  In  which  the  store  Is  Immersed.  To  reduce  the  sensitivity  to  shock  waves, 
the  pressures  at  the  top  and  bottom  of  each  section  of  the  body  are  assigned  the  values 
at  the  intersection  of  Mach  lines  as  Illustrated  in  Fig.  5.66. 


TOP  POINT 


FIG.  5.64  BUOYANCY  CORRECTION  DIAGRAM 


Then  the  vertical  buoyant  force  coefficient  per  unit  length  i»  given  by 


Cm  * - — f. 
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(5.167) 


The  pressure  Is  assumed  to  vary  sinusoidally  from  a maximum  at  the  bottom  point  to  a 
minimum  at  the  top  point.  Then  the  vertical  component  of  the  force  coefficient  results 
in  the  buoyant  force  coefficient  shown  in  Eq.  (5.167).  R is  the  local  body  radius,  and 
Sref  is  a reference  area,  ACp  Is  the  difference  In  pressure  coefficient  between  the  two 
points  and  Is  found  from  the  linear  theory  expression 


CP  * - & 


(5.168) 


where  u Is  the  local  perturbation  velocity  of  the  stream  in  which  the  store  is  placed. 

Side  forces  due  to  buoyancy  are  found  in  an  analogous  manner  using  the  lateral 
pressure  difference  across  the  store  to  obtain  ACp. 

The  buoyant  forces  are  generally  found  to  be  small  when  calculated  in  this 
manner,  and  any  moments  resulting  from  buoyancy  effects  are  neglected. 

5.6.2  Loading  on  low  aspect  ratio  store  fins 

High  aspect  ratio  fins  are  treated  differently,  while  forces  on  low  aspect 
ratio  fins  are  obtained  by  a slender  body  theory. 

Assuming  that  an  axial  section  of  the  store  contains  part  of  a low  aspect 
ratio  fin,  then  the  components  of  Eq.  (5.163)  become 
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(5.169) 


CN  “ CN  I body  + I fin  + carry-over 


(5.170) 


From  the  theory  for  slender  wing-body  interference,  given  in  ref.  84 , the  fin  contribu- 
tion (lift  panels)  is 


dCN,  8 
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(b22  k„2  - bj2  K^) 


(5.171) 


Here  bj  and  b2  are  the  fin  semi-span  dimensions  illustrated  in  Fig.  5. 45 . and 

are  wing  carry-over  factors  at  stations  x^  and  x2.  Values  are  given  in  ref.  8 4 . D is 
the  body  reference  diameter  used  in  the  non-dimensionallzation  of  the  forces. 

Also,  based  on  the  slender  body  theory  of  ref.  8^ 

Sr' carry-over  * (b2?KB2  ' bl2*V  (5‘172) 

where  Kg^  and  Kg  are  body  carry-over  factors,  al30  given  in  ref.  8^. 


FIG.  5.A5  PORTION  OF  BODY  WITH  LOW  ASPECT  RATIO  FINS 
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The  same  expression  for  slender  wind  lift  gives,  for  two  fins, 


fin  + body 
carry-over 


5?< h22K«2  + b!2KWl) 


(5.173) 


The  fin  is  truncated  after  the  axial  station  of  maximum  span. 

5.6.3  Loading  on  high  aspect  ratio  store  fins 

For  a high  aspect  ratio  wing,  the  force  on  two  wing  panels  (body  removed)  Is 
obtained  by  linear  theory.  The  coefficient  per  unit  angle  of  attack  Is  designated  CN  . 
The  Incremental  loading  on  the  body  due  to  wing  carry-over  Is  0 


dx  1 carry-over 


KBcNa 

(X2-X!) 


and 


(5.174) 


c I 

N'fin  + carry-over 


Cm  Is  determined  by  linear  theory  as  In  the  Mach  box  method,  for  example.  By  a similar 



process  JftaJ  Is  calculated  by  linear  theory  as  the  normal  force  coefficient  for  unit 
dx  da 

value  of  (37).  That  Is,  a varies  linearly  across  the  chord  of  the  wing,  being  zero 


KwCN 


<*2-*l>  & 


(5.175) 


at  the  center  of  the  interval  (X2-Xj). 

In  Fernandes'  procedure,  the  forces  on  the  store  due  to  the  interference  flow 
field  are  added  to  the  forces  due  to  store  angle  of  attack  and  motion.  This  linear 
superposition  works  well  unless  the  combined  angles  are  large. 


5.7  Trajectory  Calculation  In  Supersonic  Flow 


The  various  steps  can  be  carried  out  by  the  computer  codes  developed  In  ref3, 
62,  63,  ana  64.  The  flow  field  and  force  calculations  can  then  be  combined  with  a 
trajectory  program  to  provide  the  store  motion.  The  trajectory  program  part  of  the 
store  separation  code  of  ref.  56  might  be  adapted  to  this  purpose. 

5.8  Transonic  Flow  Field 


Transonic  flow  field  calculations  are  not  yet  as  fully  developed  as  the  sub- 
sonic and  supersonic  cases.  Two  methods  could  be  employed  to  approximate  the  transonic 
flow.  One  approach  Is  based  on  slender  body  theory;  while  the  other  makes  use  of 
transonic  similarity  laws. 


The  slender  body  theory  philosophy  is  based  on  the  linearized  expression  for 
velocity  potential 


(1-0 


♦xx  + ♦yy  + ♦zb 


0 


(5.176) 


For  a slender  body  $xx  is  small  compared  with  the  other  terms  and  the  flow  must  then 
satisfy  the  Laplace  equation  at  each  cross  section.  Since  the  Mach  number  is  no  longer 
a parameter,  formulas  based  on  slender  body  theory  can  be  applied  at  transonic  speeds. 

The  flow  around  wings  and  pylons,  however,  does  not  fit  into  this  category. 

Transonic  theories  are  primarily  directed  toward  the  prediction  of  forces  In 
uniform  flows.  Transonic  similarity  laws  can  be  applied,  however,  to  adapt  subsonic  or 
supersonic  theories  to  the  transonic  problem. 

One  possible  approach  would  use  slender  body  theory  to  obtain  the  flow  field 
about  bodies  and  similarity  law  to  find  the  flow  field  about  wings  and  pylons.  Then 
the  forces  on  the  store  could  be  determined  from  slender  body  theory;  or,  if  the  store 
is  not  sufficiently  slender,  from  transonic  similarity. 

From  the  transonic  similarity  laws,  which  may  be  found  in  ref.  85,  for  example, 
a given  configuration  is  transformed  according  to  the  relations 


B'AR'*  BAR 


(5.177) 
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Here  ( )'  denotes  the  transformed  variable. 


(5.178a) 


$ - /\  1 - 

AR  ■ aspect  ratio 

t * angle  of  attack  or  thickness  ratio 

r ■ (r+DwS 

y » ratio  of  specific  heats  * 1.4  for  air 

The  upwash  and  sidewash  angles  in  the  transformed  coordinates  will  scale  with 
t . That  is 


wj_  m w 

U'  t u 

and 

v'  . t ' v 
U'  T u 


(5.178b) 


(5.179) 


In  order  to  determine  a transonic  flow  field,  a supersonic  or  subsonic 
solution  is  first  obtained  then  transformed  to  the  desired  transonic  flow.  For  example, 
suppose  the  flow  field  is  desired  at  Mi  ■ .9  about  a wing  at  angle  of  attack  a . Then 
the  solution  is  obtained  for  a subsonic  case,  say  M„  » .7,  below  the  critical  speed. 

The  methods  of  Section  5-1  can  be  used  for  that  purpose.  Then  B'  ■ /7T9'  and  B - /75T. 
Also  r'  - 1.944  and  r - 1.176.  Hence,  from  Eq.  (5-177)  AR  - AR'  x |1  ■ .61  AH'  and 

from  Eq.  (5.178) 
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Then  the  subsonic  solution  for  a configuration  whose  aspect  ratio  is  .61 
times  the  desired  transonic  configuration  and  whose  thickness  and  angle  of  attack  are  .37 
times  the  corresponding  transonic  case  gives  the  required  similar  flow.  The  flow  angle 
at  corresponding  points  scale  with  I_  and  the  points  themselves  scale  with  B.  That  is 

T 


il  . fi_  i 

x'  B'  x 


(5.180) 


and 


z1  . 8 z 
x'  B ' x' 

5.9  Forces  on  Store  in  Transonic  Flow 


(5.181) 


The  forces  on  the  store  in  the  transonic  flow  field  can  be  calculated  by  slender 
body  theory  following  procedures  similar  to  those  in  Section  5.6.1  and  5.6.2.  Transonic 
similarity  laws  can  be  applied  also.  From  ref.  85 


, <f)2/3 

CF  ■ V 1/3  CF  C , B'AR'  ] 

(f-)1/3  ( r ’ t ' ) 1/3 

(5.182) 

(xl)2/3 

c-  \pi/3c"[7i^v5’ 

(5.183) 

where  Cp  and  Cp  represent  force  coefficients,  and  Cm  and  represent  moment  coefficients. 
5.10  Control  of  Trajectories 

Accurate  procedures  for  calculating  store  trajectories  can  be  used  to  investi- 
gate the  safety  of  aircraft  clearance  and  to  provide  predictions  of  Impact  points  for 
fire  control  information.  However,  another  means  of  achieving  these  ends  would  be  to 
adjust  the  ejection  process  so  as  to  keep  the  store  on  a ballistic  trajectory. 
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In  the  pitch  plane,  the  trajectory  might  be  controlled  by  ualng  a two-piston 
ejector,  with  the  ejection  forcee  adjusted  to  Impart  Just  the  right  amount  of  Initial 
linear  and  angular  velocity  to  the  store.  Then  the  store  would  emerge  from  the  aircraft 
flow  field  nearly  on  a ballistic  trajectory;  that  Is,  with  very  small  angle  of  attack 
and  pitch  rate.  The  store  would  then  continue  on  a predictable  ballistic  trajectory 
to  impact. 


Such  a controlled  launch  would  have  several  advantages  over  an  ejection  with 
larger  pitch  amplitudes.  For  one  thing,  the  smooth  trajectory  would  almost  certainly 
result  in  an  Improved  safety  margin  with  respect  to  avoidance  of  store-aircraft  collision. 
It  would  also  broaden  the  envelope  of  safe  launch  conditions  to  nearly  fill  the  entire 
range  of  aircraft  flight.  Furthermore,  impact  points  would  be  relatively  insensitive 
to  variations  in  ejector  cartridge  Impulse  or  departures  of  other  ejection  parameters 
from  their  nominal  values.  The  low  sei  iltivity  would  result  from  the  fact  that  several 
degrees  of  maximum  store  angle  of  attack  would  have  little  effect  on  the  trajectory  of 
a stable  store.  A controlled  ejection  could  also  keep  store  angles  of  attack  low  enough 
to  avoid  nonlinear  Instability  ‘s , such  as  catastrophic  yaw. 

The  investigations  reported  in  ref.  17  indicate  that  even  in  severe  cases  it 
is  possible  to  achieve  initial  angular  rates  that  would  keep  the  store  on  a ballistic 
trajectory.  In  fact,  by  applying  initial  yawing  moments  as  well  as  pitching  moments, 
the  trajectory  can  be  corrected  in  both  planes  as  might  be  required  for  launch  from  a 
wing  station  into  a lateral  flow  field.  A complete  description  of  the  experiments  can 
be  found  in  ref.  86. 

In  general,  the  distribution  of  ejection  forces  required  to  Insert  the  store 
on  a ballistic  trajectory  is  a function  of  the  store-aircraft  geometry  and  the  flight 
conditions.  Consequently,  the  ejection  mechanism  should  be  adjustable  in  flight  to 
produce  the  desired  store  Insertion  characteristics.  Several  methods  have  been  proposed 
for  mechanizing  the  adjustment. 

In  ref.  17  an  adjustable  ejection  system  is  pictured  with  four  pistons  capable 
of  applying  both  pitch  and  yaw  corrections.  The  adjustment  is  accomplished  by  controlling 
the  time  over  which  high  pressure  is  applied  to  each  piston.  A relief  valve  on  each 
piston  can  be  opened  to  cut  the  impulse  at  that  particular  ejection  foot.  The  timing 
of  the  valves,  in  turn,  would  be  set  from  the  fire  control  computer  which  would  take 
account  of  the  store  location  and  flight  conditions. 

Several  mechanisms  have  been  constructed  with  two-piston  ejectors  in  which 
the  relative  distribution  of  gas  pressure  to  the  pistons  is  controlled  by  an  adjustable 
orifice  The  orifice  turns  like  a stopcock  from  full  open  to  closed,  letting  proportionate 
amounts  of  high-pressure  gas  (from  a cartridge)  reach  the  controlled  piston.  In  order 
to  permit  in-flight  adjustment,  it  would  be  necessary  to  set  the  orifice  opening  v(  a 
an  output  from  the  fire  control  computer. 

It  is  also  possible  to  adjust  the  length  of  stroke  of  the  pistons,  and  a device 
for  accomplishing  this  has  been  constructed  at  the  Naval  Weapons  Center,  China  Lake. 

This  device  uses  auxiliary  chambers  with  a hydraulic  force  transfer  system.  The  volumes 
of  the  auxiliary  chambers  are  adjustable  so  that  portions  of  the  hydraulic  fluid  can 
be  directed  to  the  piston  or  the  auxiliary  chamber.  If  all  of  the  fluid  is  directed 
to  the  piston,  it  13  driven  to  its  full  extension.  If  the  auxiliary  chamber  is  partially 
open,  then  some  of  the  fluid  will  not  go  to  the  ejection  piston,  and  its  stroke  will  be 
correspondingly  reduced.  As  in  the  other  systems,  an  in-flight  adjustment  of  the 
auxiliary  chamber  volume  could  be  set  by  the  fire-control  computer. 

6.  WIND  TUNNEL  TEST  TECHNIQUES 

6.1  Introduction 

Three  types  of  wind  tunnel  test  are  in  use  for  determining  store  separation 
trajectories.  In  general,  only  the  part  of  the  store's  flight  in  the  vicinity  of  the 
aircraft  is  investigated.  Once  the  store  leaves  the  effective  region  of  Influence  of 
the  launch  aircraft,  it  is  simply  an  aerodynamic  shape  in  a uniform  flow  whose  motion 
can  be  determined  by  conventional  means  if  the  angle  of  attack  is  not  excessive. 

One  type  of  test,  obviously,  is  to  perform  a simulated  drop  in  the  wind  tunnel. 
The  scale  is  usually  rather  small,  and  compressibility  effects  cannot  be  correctly 
simulated  because  gravitational  acceleration  is  fixed  and  dictates  the  time  scale  once 
the  length  scale  ha3  been  chosen. 

Another  test  method  is  the  "captive  trajectory"  technique  in  which  the  store 
model  is  mounted  on  a balance  and  placed  in  its  launch  position  near  a model  of  the 
launch  aircraft.  The  forces  and  moments  or.  the  store  are  measured.  Then  a computer 
calculates  the  motion  of  the  store  during  a short  time  interval  taking  account  of  the 
measured  aerodynamic  loads,  the  gravitational  force,  and  damping  forces.  The  store  is 
then  moved  to  the  next  predicted  location  where  the  process  is  repeated.  Modern  computers 
are  fast  enough  to  carry  out  a complete  trajectory  in  a few  seconds. 

A third  commonly  used  technique  is  the  "grid  data  bank"  procedure.  Here  the 
aircraft  flow  field  is  surveyed  using  the  store  model  as  a probe,  at  every  location  the 
forces  on  the  store  due  to  its  angle  of  attack  and  the  aircraft  flow  field  are  measured 
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and  stored  In  a computer  memory.  Then  a trajectory  can  be  computed  for  arbitrary  Initial 
conditions  using  Interpolations  of  the  data  to  determine  the  forces  at  arbitrary  loca- 
tions. Gravitational  and  damping  forces  must  be  added  In  the  trajectory  program.  In 
most  cases,  the  forces  due  to  aircraft  flo*  field  and  store  angle  of  attack  can  be 
linearly  superimposed  so  that  during  the  flow  survey  the  store  can  be  maintained  at  zero 
angle  of  attack.  The  store  static  and  dynamic  characteristics  then  are  measured  only 
once  In  a uniform  stream. 

A variation  of  the  captive  trajectory  technique  uses  a flow  angularity  probe 
rather  than  a store  model.  The  probe  measures  the  flow  field  In  the  region  that  would 
be  occupied  by  the  store  axis.  Based  on  the  survey  data,  a computer  calculates  the 
forces  on  a store  and  predicts  an  Increment  of  Its  motion.  The  probe  Is  then  moved  to 
the  new  location  where  the  process  Is  repeated  until  the  trajectory  Is  traced  out. 
Although  cumbersome  In  some  ways,  this  technique  has  the  advantage  of  requiring  only 
that  a smail  probe  be  placed  In  the  flow  field,  so  that  Interference  between  the  probe 
support  and  the  aircraft  model  can  be  avoided. 

6.2  Dynamic  Drop 

6.2.1  Scaling  laws 

6. 2. 1.1  Scaling  problem 

A simple  direct  technique  for  obtaining  trajectory  data  Is  the  ejection  of  a 
scale  model  In  a wind  tunnel.  Special  equipment  Is  required  to  eject  the  model  In 
addition  to  high-speed  photography,  preferably  capable  of  recording  both  pitch  and  yaw 
motions.  Descriptions  of  tests  of  this  type  may  be  found  in  refs.  87-91.  Since  an 
airplane  model  Is  required  In  addition  to  the  store,  the  scale  of  the  test  Is  usually 
quite  small.  Consequently,  scaling  laws  are  of  some  significance  and  will  be  briefly 
reviewed  here. 

Assume  that  subscript  ( )f  refers  to  the  full-Bcale  condition  and  ( )m  denotes 
the  wind  tunnel  model.  The  objective  Is  to  obtain  a trajectory  in  the  wind  tunnel  which 
Is  In  exact  proportion  to  the  full-scale  case.  That  is,  the  store  model  will  trace  the 
same  path  with  respect  to  the  airplane  and  will  assume  identical  angular  orientations. 

6. 2. 1.2  Length 

Define  the  length  scale,  then,  by 


Lf  “ *«* 


(6.1) 


6. 2. 1.3  Time  and  velocity 

In  the  absence  of  aerodynamic  forces,  the  store  will  accelerate  at  the 
gravitational  rate.  To  make  the  displacement  of  the  model  scale  In  the  same  manner  as 
Its  own  length,  the  tine  will  scale  according  to 
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and  the  velocity  scale  will  be 
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6. 2.1. A Air  temperature 

In  order  to  correctly  simulate  compressibility  effects,  the  Mach  number  '"■ust  be 

matched 


(6.14) 


Hence  the  velocity  of  sound  Is  constrained  as  follows 
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Since,  for  a perfect  gas  the  velocity  of  sound  Is  proportional  to  the  square  root  of  the 
absolute  temperature 
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In  most  supersonic  wind  tunnels,  this  requirement  cannot  be  realized.  The 
model  is  so  small  that  the  air  would  condense.  The  consequences  will  be  discussed  later. 

6. 2. 1.5  Weight  scale 

To  make  aerodynamic  forces,  P,  scale  in  proportion  to  force  of  gravity, 


m 

Pf 
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The  aerodynamic  forces,  in  turn,  are  proportional  to  vehicle  area  and  stream  dynamic 
pressure.  Hence,  for  a properly  scaled  test 
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m _ 1 Km 
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where  pf  is  the  density  of  the  air  in  the  free  flight  environment,  and  pm  is  the  density 
of  the  air  in  the  wind  tunnel  teat  section. 


If  Mach  number  is  matched,  even  if  the  wind  tunnel  velocity  is  not  properly 
scaled,  then  the  required  weight  scale  becomes 


— . 1 
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where  pf  and  pm  are  the  static  pressures  in  free  flight  and  the  wind  tunnel,  respectively. 
This  result  is  a consequence  of  the  fact  that  dynamic  pressure  is  related  to  static 
pressure  in  compressible  flow  of  a perfect  gas  by 
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Y being  the  ratio  of  specific  heats. 

6.2. 1.6  Moment  of  Inertia 


In  order  to  make  the  model  turn  through  the  same  angle  at  the  same  point  on 
the  trajectory  as  the  full-scale  vehicle,  the  moments  of  inertia  must  be  related  by 
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By  comparison  with  the  weight  scale  from  Eq,  (6.9)  the  radii  of  gyration  will  have  the 
ratio 
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6. 2. 1.7  Aerodynamic  damping 

In  order  to  properly  scale  the  aerodynamic  damping,  the  velocity  of  the  stream 
should  scale  in  the  3ame  ratio  as  the  trajectory  velocities  as  given  by  Eq.  (6.3).  If 
the  Mach  numbers  are  matched,  the  wind  tunnel  stream  will  have  too  high  a velocity,  Ura, 
and  hence  the  model  store  will  not  have  the  correct  damping. 

6. 2. 1.8  Induced  angle  of  attack 

The  supersonic  wind  tunnel  with  its  incorrect  velocity  will  also  have  the  wrong 
Induced  angle  of  attack. 

6. 2. 1.9  Choice  of  scaling  laws 

Unless  the  tunnel  is  so  large  that  nearly  full-scale  models  can  be  tested,  it 
is  not  possible  to  match  all  non-dimensional  parameters.  It  is  possible,  however,  to 
choose  which  parameter  will  not  be  properly  scaled.  The  following  table  shows  turee  sets 
of  scaling  laws, 
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The  heavy-model  scaling  of  set  I produces  a trajectory  nearly  geometrically 
similar  to  the  free-fllght  trajectory  and  Is  generally  preferred  in  supersonic  tests 
unless  the  store  acquires  a large  vertical  velocity.  In  that  case,  the  Induced  angle 
of  attack  will  be  too  small  In  the  model  test  and  a noticeable,  and  usually  uncor.serva- 
tlve,  discrepancy  In  the  trajectory  may  result. 

Set  II  Is  satisfactory  at  low  speeds  when  Mach  number  Is  not  a significant 
parameter.  It  is  not  reliable  at  supersonic  speeds,  however. 

In  set  III  the  aerodynamic  accelerations  acting  on  the  model  would  be  different 
from  the  full  scale,  but  gravity  will  not  oblige  by  changing  In  the  same  ratio. 
Nevertheless,  this  set  of  scaling  laws  is  appropriate  for  cases  In  which  a substantial 
vertical  ejection  velocity  Is  Initially  applied  to  the  store  so  that  the  displacement 
due  to  gravity  Is  not  significant  In  the  vicinity  of  the  airplane.  Corrections  may  be 
possible.  The  measured  vertical  displacement  can  be  corrected  to  account  for  the 
discrepancy  In  gravity.  Also,  it  Is  technically  possible  to  adjust  gravity  by  applying 
a suitable  magnetic  field  (see  ref.  92).  It  is  also  conceivable  that  the  launch  air- 
plane might  be  moved  vertically  during  the  drop  test  in  order  to  correct  the  relative 
displacement . 

6.2.2  Comparison  with  flight  test 

A comparison  of  wind  tunnel  and  flight  test  releases  of  an  empty  rocket  pod 
Is  shown  in  Fig.  6.1  taken  from  ref.  93.  Although  the  camera  positions  were  different, 

It  is  apparent  that  the  wind  tunnel  test  gave  a good  representation  of  the  store 
trajectory.  Comparisons  of  vertical  position  and  pitch  angle  for  a f in-stabilized  bomb 
are  shown  in  Fig.  6.2  taken  from  ref.  88. 

In  the  dynamic  drop  test  technique,  one  of  the  most  difficult  problems  is  the 
calibration  of  the  ejector  system  to  cause  it  to  provide  the  desired  initial  conditions. 
Eventually  it  must  be  calibrated  by  direct  measurement  of  the  initial  linear  and  angular 
velocity  imparted  to  stores  under  actual  wind-on  conditions. 

6.2.3  Data  reduction 

The  dynamic  drop  data  appears  in  the  form  of  photographic  coverage  of  the 
trajectory  from  high-speed  movie  film.  Reading  the  film  is  a tedious  process,  and 
accuracy  is  limited.  Procedures  have  been  developed  which  speed  the  process  and  improve 
accuracy . 

The  method  reported  In  ref.  91*  projects  the  photographic  data  on  a television 
picture  of  a model  of  the  store.  The  model  orientation  Is  adjusted  to  match  the 
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trajectory  photograph.  Then  the  settings  that  produce  the  match  are  automatically 
recorded.  This  process  was  developed  primarily  to  reduce  flight  test  film  data  of 
store  separation  experiments. 

6.3  Captive  Trajectory 

In  the  captive  trajectory  technique,  an  airplane  model  Is  placed  In  the  wind 
tunnel,  and  a store  model  is  mounted  on  a separate  sting  balance  which  places  it  close 
to  Its  carriage  position.  The  forces  on  the  store  are  measured  by  Its  balance.  These 
forces  are  fed  Into  a computer  along  with  Inertial  and  damping  data.  The  computer 
calculates  an  Increment  of  motion  which  Is  then  fed  to  the  support  mechanism  which 
automatically  moves  the  store  to  the  next  position  along  Its  trajectory.  The  result  Is 
a "captive  trajectory"  In  which  the  model  store  follows  the  same  path  as  a full-scale 
store  would  travel  under  the  same  initial  conditions.  Such  systems  are  described  in 
refs.  95-9/. 

The  technique  seems  to  work  well.  Comparison  with  flight  test  results  Is 
difficult  because  flight  conditions  and  wind  tunnel  conditions  are  seldom  matched.  How- 
ever, wind  tunnel  trajectories  agree  with  flight  test  when  valid  comparison  Is  possible. 
Figures  6.3  and  6. A from  ref.  98,  for  example,  compare  captive  trajectory  and  flight 
test  data  for  a Jettison  of  a gun-pod  from  an  F-A  aircraft. 

Difficulties  arise  in  starting  the  store  very  close  to  the  carriage  position 
without  coercing  the  balance  by  contact  between  the  store  model  and  the  airplane  model. 
Also,  as  the  store  pitches  around  under  the  airplane,  it  sometimes  becomes  Impossible 
to  find  a means  of  supporting  the  store  without  having  the  support  run  into  the  air- 
plane model.  In  fact,  the  support  mechanism  causes  significant  blockage  and  flow 
interference  that  can  modify  the  trajectory.  These  and  other  factors  affecting  the 
accuracy  of  the  technique  are  discussed  In  ref.  99. 

A modification  of  the  captive  trajectory  procedure  makes  use  of  a flow-probe 
Instead  of  the  store  model.  The  probe  surveys  trie  flow  along  the  axis  of  the  store 
position.  Then  this  flow  data  is  used  to  calculate  the  force  on  the  store  In  that 
position.  Again  a trajectory  computation  Is  used  to  fix  successive  positions  of  the 
store  which  are  then  surveyed  In  turn.  The  use  of  this  technique  Is  described  In 
ref.  100. 


6 . A Grid  Data  Bank 

This  flow  survey  procedure  uses  an  airplane  model  mounted  In  a wind  tunnel  as 
in  the  captive  trajectory  technique.  However,  instead  of  moving  a store  model  along 
Its  computed  trajectory,  the  store,  mounted  on  a sting  balance,  is  used  as  a probe  to 
measure  the  aircraft  Interference  flow  field  (ref.  101).  Once  the  forces  and  moments 
on  the  store  have  been  measured  at  a grid  of  points  in  the  aircraft  flow  field,  the 
data  is  stored  In  the  form  of  Interference  influence  coefficients  in  a computer.  The 
store  static  and  dynamic  characteristics  as  a function  of  angle  of  attack  are  also 
measured  in  a uniform  stream. 

To  calculate  a trajectory,  an  initial  condition  puts  the  store  at  seme  point 
in  the  flow.  The  store-alone  forces  due  to  its  motion  and  angle  of  attack  arc  then 
added  to  the  interference  forces  interpolated  from  the  data  bank,  and  a resulting 
Increment  of  motion  is  computed.  The  computer  thus  steps  through  the  trajectory  using 
Interpolated  Interference  data  plus  store  alore  coefficients  at  each  increment. 

Once  the  flow  field  has  been  surveyed  using  the  store  for  a probe,  any  number 
of  initial  conditions  can  ue  Investigated.  However  a separate  survey  Is  required  for 
each  aircraft  angle  of  attack  ar.d  each  store  geometry. 

Comparisons  of  this  type  of  trajectory  with  dynamic  drop  data  are  shown  In 
Figs.  6. 5-6. 8,  taken  from  ref.  17.  The  store  model  represented  a 900-pound  bomb 
ejected  from  an  F-lA  aircraft  wing  pylon.  q„  and  r0  are  the  initial  angular  rates  of 
the  store  in  pitch  and  yaw,  v0  and  v0  are  the  corresponding  linear  initial  velocities. 

P y 

The  test  was  carried  out  at  Mach  number  2.5.  The  curves  marked  "experimental  data" 
represent  the  dynamic  drop  trajectories,  while  the  "calculated  data"  used  the  grid  data 
bank  technique. 

7 . FLIGHT  TEST 

7.1  Introduction 

Flight  tests  finally  provide  the  real  store  separation  characteristics  which 
analytical  and  ground  test  methods  are  trying  to  predict.  However,  while  flight  test 
data  may  be  indisputable.  Its  proper  interpretation  is  another  matter.  The  first 
problem  is  to  be  sure  of  safe  separation  before  the  test  is  undertaken. 

The  next  step  is  to  acquire  trajectory  data.  The  best  information  is  obtained 
from  high-speed  motion  pictures  from  boresighted  on-board  cameras.  Ground-based  cameras, 
or  photography  from  nearby  aircraft  can  also  be  helpful.  Converting  the  film  record 
into  a trajectory  is  a task  of  significant  dimension,  from  t;,e  points  of  view  of  both 
quality  and  quantity. 


Z DISPLACEMENT  (VERTICAL) 


FIG.  6.3  COMPARISON  WITH  FLIGHT  TEST  OF  VERTICAL  DISPLACEMENT  OF  STORE 
MEASURED  BY  CAPTIVE  TRAJECTORY  TECHNIQUE 
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FIG.  6.4  COMPARISON  WITH  FLIGHT  TEST  OF  PITCH  ANGLE  OF  STORE 
MEASURED  BY  CAPTIVE  TRAJECTORY  TECHNIQUE 


C.G.  LOCATION  HORIZONTAL 
(CALIBERS) 


FIG.  6.7  COMPARISON  OF  MEASURED  AND  CALCULATED  BOMB  MODEL 
CENTER  OF  GRAVITY  MOTION  IN  THE  YAW  PLANE 
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FIG.  6.8  COMPARISON  OF  MEASURED  AND  CALCULATED  BOMB 
MODEL  ANGLE  OF  SIDESLIP,  B,  HISTORY 
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The  most  difficult  part  of  the  flight  test  procedure,  however,  is  maintaining 
and  measuring  aircraft  flight  conditions. 

Each  of  these  problem  areas  will  be  discussed  briefly.  New  techniques  are 
constantly  being  developed,  but  implementation  is  slow  because  of  the  cost  and  hazardous 
nature  of  flight  hardware  Innovations. 

7.2  Safe  Separation 

Safe  separation  may  be  assured  by  three  methods.  First,  Covert's  criterion 
(see  Section  9)  or  other  simple  estimates  of  initial  store  motion  (such  as  experience 
with  similar  cases)  can  be  applied  to  distinguish  safe  from  hazardous  situations.  The 
advantage  of  such  procedures  lies  In  the  ability  to  rapidly  select  appropriate  release 
conditions  for  further  study.  Next,  full  trajectory  calculations  can  be  made  by  use  of 
the  more  detailed  methods  described  in  Section  5.  Finally  wind  tunnel  tests  can  be 
made  and  compared  with  the  analytical  predictions. 

The  first  flight  tests  of  a new  store-aircraft  combination  naturally  would  be 
made  at  a very  safe  condition.  That  is  one  in  which  wind  tunnel  tests  and  analytical 
predictions  are  in  agreement,  and  for  which  theory  and  experiment  Indicate  a wide  margin 
of  safety  in  case  flight  or  release  conditions  depart  somewhat  from  the  desired  values. 

Once  flight  test  data  becomes  available,  it  can  be  compared  to  analytical  and 
wind  tunnel  predictions.  Any  significant  discrepancies  must  be  resolved  before  flight 
tests  are  resumed.  Assuming  that  the  flight  data  can  be  accurately  correlated  with 
ground  test  and  analysis,  the  flight  program  continues  by  safe  increments  toward  the 
boundaries  of  the  3afe-separation  or  airplane  performance  envelop.  The  number  of 
possible  store-aircraft  flight-condition  combinations  is  extremely  large;  hence  it  is 
not  feasible  to  test  every  case.  Therefore  the  extreme  conditions  are  identified  and 
tested,  thus  bounding  a domain  of  safe  separation. 

7.3  Instrumentation 

For  proper  Interpretation  of  store  separation  data,  it  is  necessary,  of 
course,  to  know  the  flight  condition  of  the  launch  aircraft.  This  information  is 
obtained  by  the  flight  instrumentation  conventionally  available,  but  unusual  accuracy  is 
required  in  some  measurements.  The  information  required  to  establish  the  flight 
condition  of  a particular  test  flight  Includes 

flight  speed  and  Mach  number 
altitude  (air  density  and  temperature) 
dive  angle 

roll  and  yaw  angles,  if  any 
angle  of  attack 
maneuver,  if  any 

In  cases  of  steady  flight  in  a vertical  plane  (no  yaw  or  roll),  the  various 
flight  parameters  can  be  obtained  with  conventional  aircraft  instruments.  To  assure 
steady  flight,  accelerometers  can  monitor  inadvertent  maneuvers. 

Since  store  separation  trajectories  are  not  critically  sensitive  to  Mach 
number,  dive  angle,  or  altitude,  measurement  of  these  quantities  is  not  difficult. 

Sudden  aircraft  maneuvers  due  to  turbulence,  for  example,  could  disturb  the  store 
trajectory  or  impart  motions  to  cameras  photographing  the  store  separation  which  might 
be  interpreted  as  store  motion.  Consequently,  for  determination  of  ballistics  or  basic 
separation  characteristics,  steady  flight  is  essential. 

The  airplane  angle  of  attack  must  be  accurately  known  to  fully  define  the 
separation  data.  It  can  be  measured  by  a flow  probe,  but  probably  the  most  accurate 
determination  would  use  known  aircraft  weight,  flight  speed,  dive  angle,  altitude,  and 
lift  coefficient  slope  to  calculate  angle  of  attack  (during  zero  acceleration). 

Reference  102  briefly  describes  the  use  and  calibration  of  aircraft  flight 
instrumentation.  We  are  concerned  here  primarily  with  the  problem  of  measuring  the 
trajectory  of  the  store  with  respect  to  the  airplane.  The  initial  condition  imparted 
by  the  ejection  system  is  required,  as  well  as  the  resulting  trajectory  of  the  store. 
Such  information  as  ln-carriage  loads,  structural  and  aeroelastic  effects,  and  effects 
of  the  store  on  airplane  performance  may  also  be  of  interest,  but  are  not  essential  to 
the  determination  of  the  trajectory  which  is  the  item  under  consideration. 

Exhaustive  tests  of  ejector  characteristics  are  described  in  ref.  10’  They 
indicate  surprisingly  reliable  and  repeatable  performance  of  cartridge-powered  ejectors. 
However,  because  of  occasional  vagaries,  detailed  examination  of  store  trajectories 
should  make  use  of  the  observed  Initial  linear  and  angular  rates  imparted  to  the  store 
as  well  as  the  calibration  of  the  ejection  system. 

Generally,  then,  the  trajectory  data  of  Interest  is  provided  by  photographic 
coverage  from  three  sources.  Ground-based  theodolites  provide  accurate  measurement  of 
store  and  aircraft  position  as  functions  of  time  from  well-established  ground  reference 
points.  Unfortunately,  the  distance  is  too  great  to  give  good  data  on  store  orientation 
and  position  with  respect  to  the  airplane. 


97 


Photographic  Information  from  chaBe  planes  Is  of  more  quantitative  value 
because  the  airplane  can  get  close  enough  to  the  separation  process  to  provide  good  data 
on  the  motion  of  the  store  relative  to  the  launch  aircraft.  The  chase  plane  camera, 
however.  Is  not  accurately  located  with  respect  to  the  launch  airplane,  and  hence  the 
data  reduction  process  Is  only  as  accurate  as  the  estimate  of  chase  plane  position. 

The  chase  plane  can  al3o  be  located  by  ground  theodolites,  and  could  even  be  photo- 
graphed from  fixed  cameras  on  the  launch  airplane  to  triangulate  its  relative  position. 
However,  the  ground  coverage  again  suffers  because  of  Its  distance  from  the  event. 
Triangulation  of  the  chase  plane  from  the  launch  aircraft  Is  also  Impractical  unless  the 
chase  airplane  lies  within  the  predetermined  field  of  view  of  triangulating  cameras. 

Consequently,  the  most  accurate  data  on  the  separation  process  and  the  Initial 
store  trajectory  Is  obtained  from  carefullv  boreslghted  cameras  mounted  on  the  launch 
airplane.  The  following  accuracies  are  desired,  but  may  not  be  required  In  all  cases 

a)  position  of  points  on  the  3tore  to  .1  foot  with  respect  to  aircraft 
coordinates 

b)  pitch,  yaw,  and  roll  attitude  to  1 degree  relative  to  aircraft  coordinates 

c)  store  control  deflections  to  1 degree 

d)  time  reference  data  accurate  to  .001  second 

A portion  of  the  aircraft  fuselage  should  be  In  the  field  of  view  to  serve  as 
a reference  for  aircraft  structural  deformation  corrections. 

Ejector-foot  time  history  measurement  In  flight  would  also  help  to  establish 
Initial  conditions  provided  by  the  ejection  process. 

Accuracies  actually  achievable  by  a photographic  system  have  been  reported  In 
ref.  109.  In  flight,  the  store,  marked  with  special  reference  Indicators,  is  released 
and  photographed  by  on-board  cameras.  The  film  is  later  processed,  frame  by  frame,  and 
the  locations  of  the  reference  marks  are  stored  on  computer  cards.  The  data  cards 
plus  the  camera  and  reference  mark  locations  are  fed  into  a computer  which  i3  programmed 
to  resolve  the  geometric  relationships  and  calculate  store  trajectory  information. 

To  test  the  accuracy  of  the  procedure,  cameras  were  mounted  on  fixed  platforms 
(on  the  ground)  and  a store  wa.,  placed  at  various  carefully  measured  locations. 
Photographic  data  was  then  processed  in  the  same  manner  as  flight  test  data,  and  the 
calculated  store  positions  were  compared  with  the  measured  ones.  Samples  of  the 
resulting  discrepancies  are  shown  in  Figs.  7.1  and  7.2  taken  from  ref.  109.  From  these 
results,  it  Is  apparent  that  the  angular  errors  are  more  like  2 degrees  than  the  desired 
1 degree,  and  displacement  error  Is  also  several  times  as  large  as  the  desired  .1  foot. 

7. A Data  Reduction 

Knowing  the  positions  and  optical  axes  of  two  aircraft-mounted  cameras,  It  is 
possible  to  determine  the  precise  location  of  points  on  a separating  store  seen  by  both 
cameras.  Reducing  the  data  normally  requires  tedious  measurements  of  coordinates  of  the 
desired  points  from  fiducial  marks  on  frames  of  movie  film.  The  data  can  be  reduced  by 
a computer  once  the  positions  have  been  recorded;  however  several  sources  of  error 
limit  the  accuracy  of  the  process.  The  resolution  and  repeatability  of  the  film  reading 
process  provides  a fundamental  limitation  on  the  accuracy.  In  addition,  distortions 
produced  by  camera  lenses  and  films  can  cause  errors  although  corrections  are  possible. 
Another  source  of  error  Is  the  displacement  in  time  between  the  corresponding  movie 
frames  taken  by  two  different  cameras.  As  indicated  by  the  experiments  of  ref.  109, 
these  errors  are  likely  to  produce  Inaccuracy  In  angle  measurements  of  about  2 degrees, 
and  displacement  discrepancies  of  several  Inches.  Structural  deformation  of  the  air- 
craft can  cause  camera  motions  that  will  Introduce  additional  errors. 

Several  procedures  have  been  investigated  to  reduce  the  labor  and  improve 
the  accuracy  of  the  data  reduction  process.  A method  called  "Teledaq"  (Television  Data 
Acquisition)  Is  proposed  in  ref.  105.  In  this  scheme,  television  cameras  are  used 
instead  of  movie  cameras,  so  that  the  data  Is  automatically  available  in  the  form  of 
electrical  signals  on  tape.  To  avoid  the  difficulties  of  pattern  recognition,  the 
store  is  to  be  painted  red  In  certain  areas,  and  a red  filter  on  part  of  the  camera 
system  is  used  to  make  the  red-painted  areas  particularly  dominant.  With  this  type  of 
coding,  the  television  tapes  can  be  digitized  and  the  store  position  and  orientation 
determined  by  computerized  operations  on  the  data.  Such  a system  could  greatly 
facilitate  data  reduction.  According  to  estimates  given  in  ref.  105,  the  accuracy  would 
be  about  the  same  as  photographic  methods,  but  advances  being  developed  in  data  storage 
and  resolution  could  lead  to  significantly  Improved  versions. 

Another  procedure,  described  previously  as  being  applicable  to  the  reduction 
of  wind  tunnel  data,  has  been  developed  at  the  Naval  Missile  Center  (ref.  99).  The 
FDAG  (Photo  Data  Analysis  System)  uses  aircraft-carried  cameras  to  record  the  fllyht, 
but  reads  the  film  by  superimposing  the  image  of  a store  model  viewed  from  a television 
camera  and  an  aircraft  camera  frame  showing  the  actual  store.  When  the  model  has  been 
correctly  positioned  with  respect  to  the  television  camera  so  that  it  exactly  coincides 
with  the  movie  camera  image,  then  the  model  position  and  orientation  is  automatically 
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recorded.  The  trajectory  data  can  be  accurately  obtained  In  this  manner  with  only  one 
on-board  camera.  However,  accuracy  Is  Improved  by  statistically  averaging  data  from 
two  or  more  cameras. 

By  smoothing  the  data  obtained  by  the  PDAS  system,  It  Is  possible  to  differ- 
entiate the  position  Information  to  obtain  velocities  and  accelerations  of  the  store. 
Consequently  the  forces  on  the  store  can  be  deduced  for  comparison  with  analytical 
predictions  or  wind  tunnel  measurements. 

A procedure  to  Improve  angular  resolution  of  the  store  is  proposed  in  ref. 

106.  The  Idea  is  to  place  carefully  oriented  polarized  reflectors  on  the  store.  When 
these  are  Illuminated  with  polarized  light,  the  angle  at  which  the  reflected  light  is 
extinguished  can  be  determined,  and  hence  the  orientation  of  the  store  at  that  instant 
is  known.  Other  possible  uses  of  polarized  light  sources  are  also  Indicated  In  the 
report  (ref.  106). 

8.  CONCLUDING  REMARKS 

The  primary  subject  of  this  report  has  been  the  trajectories  of  external 
stores  during  their  passage  through  the  flow  field  of  the  carrying  aircraft.  Easily 
applied  safe-separation  criteria  are  described  In  Section  A,  and  methods  of  calculating 
trajectories  are  indicated  In  Section  5. 

Detailed  numerical  calculations  of  subsonic  trajectories  can  be  carried  out 
using  the  computer  programs  described  In  refs.  56-58.  he  analysis  seems  to  be  quite 
accurate  and  Is  capable  of  handling  very  general  configurations  and  flight  conditions. 

The  main  limitations  are  In  angle  of  attack  and  Mach  number.  If  the  store  angle  is  such 
that  nonlinearities  due  to  body  vortex  separation  become  significant,  then  the  approxi- 
mate treatment  given  In  refs.  56-58  may  be  inadequate.  Also,  the  compressibility 
correction  procedure  becomes  unworkable  at  speeds  above  the  critical  Mach  number. 

The  supersonic  case  Is  analyzed  satisfactorily  by  Fernandes  In  refs.  62-6A. 
However,  he  does  not  include  the  trajectory  calculations  in  his  computer  programs.  His 
procedure  may  also  be  inadequate  for  stores  at  high  angles  of  attack  where  vortex 
separation  becomes  Important.  At  low  supersonic  Mach  numbers,  the  procedure  will  break 
down  because  of  the  inapplicability  of  the  linearization  of  the  flow  equations. 

Thus  the  transonic  regime  Is  the  area  that  has  been  lea3t  developed.  Some 
procedures  are  suggested  In  Sections  5.8  and  5.9,  but  comprehensive  computer  codes 
covering  the  transonic  regime  do  not  yet  appear  to  be  available. 

Although  the  effect  of  external  stores  on  airplane  performance  has  been  only 
briefly  discussed  In  this  report,  the  Implication  seems  clear  that  flight  characteristics 
could  be  significantly  Improved  by  better  design  and  Integration,  as  in  the  conformal 
carriage  concept. 

Supersonic  release  of  stores  has  received  relatively  little  attention  to  date. 
Hot  many  aircraft  can  carry  external  stores  at  supersonic  speeds,  and  the  large  aero- 
dynamic forces  make  store  trajectories  difficult  to  control.  However,  Improved  aircraft 
performance  and  the  need  for  high-speed  attacks  are  likely  to  result  in  increasing 
release  velocities.  To  assure  aircraft  safety  and  accurate  placement  of  store  trajecto- 
ries, better  ejection  mechanisms  will  be  needed  at  the  higher  speeds.  Controlled 
releases,  such  as  those  described  In  Section  5.10,  will  probably  be  developed. 

The  major  message  of  the  considerable  work  in  the  field  of  store  separation  is 
that  better  Integration  of  aircraft  and  stores  is  needed  to  improve  the  performance  and 
reliability  of  both.  Aircraft  should  be  designed  to  carry  stores,  instead  of  adding 
this  capability  almost  a3  an  afterthought.  It  should  also  be  noted  that  If  a new  air- 
plane must  be  capable  of  carrying  all  existing  stores,  and  a new  store  must  fit  on  all 
existing  aircraft,  Improvements  In  performance  of  the  aircraft-store  combination  will 
be  extremely  slow.  Integrated  store-aircraft  systems  must  be  conceived  without  such 
constraints  to  provide  paths  for  far-reaching  development. 
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